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ABSTRACT 

A spectral survey of IRC+10216 has been carried out in the range 11 to 14 /zm with a spectral 
resolution of about 4 km s _1 . We have identified a forest of lines in six bands of C2H2 involving the 
vibrational states from the ground to 2>v§ and in two bands of HCN, involving the vibrational states 
from the ground up to 2v2- Some of these transitions are observed also in H 13 CCH and H 13 CN. We 
have estimated the kinetic, vibrational, and rotational temperatures, and the abundances and column 
densities of C2H2 and HCN between 1 and 300 R» (~ 1.5 x 10 16 cm) by fitting about 300 of these 
ro-vibrational lines. The envelope can be divided into three regions with approximate boundaries at 
0'.'019 (the stellar photosphere), O'.'l (the inner dust formation zone), and 0'.'4 (outer dust formation 
zone). Most of the lines might require a large microturbulencc broadening. The derived abundances 
of C2H2 and HCN increase by factors of 10 and 4, respectively, from the innermost envelope outwards. 
The derived column densities for both C2H2 and HCN are ~ 1.6 x 10 19 cm" 2 . Vibrational states up to 
3000 K above ground are populated, suggesting pumping by near-infrared radiation from the star and 
innermost envelope. Low rotational levels can be considered under LTE while those with J > 20 — 30 
are not thermalized. A few lines require special analysis to deal with effects like overlap with lines of 
other molecules. 

Subject headings: line: identification — line: profiles — surveys — stars: AGB and post-AGB - 
stars: carbon — stars: individual (IRC+10216) 



1. INTRODUCTION 

IRC+10216 is a carbon AGB star surrounded by a cir- 
cumstellar envelope (CSE) of gas and dust. Since it is the 
nearest AGB star of this type and the strongest infrared 
object in the sky, it has become the paradigm for this 
kind of sources (see §3.1 for more information about pre- 
vious work). The central star itself cannot be observed 
directly due to the large dust optical depth but the ob- 
servational data from the CSE suggest that it is a phys- 
ically and chemically rich environment. Although this 
source has been studied many times since its discovery 
by Neugebauer & Leighton (1969), most of these physical 
and chemical processes remain poorly understood across 
the envelope (see, for example, Agundez & Cernicharo 
2006). 

The physical conditions in the innermost part of the 
CSE maintain chemical thermodynamical equilibrium. 
However, the gas could depart from the latter due to 
the action of periodic shock waves arising from stellar 
pulsation (Cherchneff et al. 1992; Agundez & Cernicharo 
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2006). At the temperatures prevailing in these zones, 
most atoms are integrated into molecular species, pri- 
marily H2. The chemistry is so rich that more than 
50% of the molecules known in space were first discov- 
ered in this source (see, for example, Morris et al. 1975; 
Betz 1981; Cernicharo et al. 2000, 2004, and references 
therein). Some of these molecules, such as SiO (Schoier 
et al. 2006), SiS, CS (Lucas et al. 1995), and metal- 
bearing species (Cernicharo & Guelin 1987; Cernicharo 
et al. 2000), are refractory and are formed in the in- 
ner envelope. As soon as the temperature of the gas 
is below a critical value these species start to condense 
and form dust grains. Other species, mainly radicals, 
display emission from the external shell of the envelope 
(Ay ~ 1 mag) (Dayal & Bieging 1993, 1995; Lucas et 
al. 1995; Cernicharo & Guelin 1996a; Lindqvist et al. 
2000) where the Galactic UV field starts to photodisso- 
ciate the stable molecules formed in the inner CSE and 
neutral-radical reactions produce long carbon chain rad- 
icals (Cernicharo & Guelin 1996a; Guelin & Cernicharo 
1997) and also the interesting gas-phase oxygen chem- 
istry (Agundez & Cernicharo 2006). 

The most abundant molecular species, after H2, are 
CO, C 2 H 2 , and HCN. CO has an abundance of ~ 8xl0" 4 
from 1 R* to beyond 1000 R„ (Lafont et al. 1982; Knapp 
& Morris 1985; Agundez & Cernicharo 2006). C 2 H 2 , 
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is predicted by thermodynamical chemical equilibrium 
models (Tejero & Cernicharo 1991) to be the most abun- 
dant molecule after CO and H 2 . Its abundance has been 
derived to be ~ 8 x 1CP 5 in the zone 1 — 40 R* from in- 
frared observations (Keady & Ridgway 1993; Cernicharo 
et al. 1999). The abundance of HCN ~ 3 - 4 x 1CT 5 has 
been derived from mid-, far-infrared, and radio observa- 
tions (Keady & Ridgway 1993; Cernicharo et al. 1996b, 
1999). 

Radiation pressure by stellar photons accelerates the 
dust grains formed near the star. For the densities pre- 
vailing in these dusty regions of the CSE, the coupling 
between gas and dust is high and hence the gas is also 
accelerated (Gilman 1972; Kwok 1975). Keady et al. 
(1988) have found that the terminal velocity of the gas 
is about 14 km s^ 1 (see also Cernicharo et al. 2000) and 
is reached before 20 R*. These authors also found that 
the velocity field seems to have more than one acceler- 
ation regime occurring at different places of the inner 
CSE. These zones might be related to the condensation 
temperature of different refractory molecular species. 

The main goal of this work is to learn more about 
the physical conditions and chemical composition of the 
CSE of IRC+10216 from modeling the line profiles of 
C2H2 and HCN observed in the mid-infrared with the 
high-resolution spectrograph TEXES (Lacy et al. 2002). 
These observations are presented in §2. The models we 
have used to fit the observational data are described in §3 
and Appendices A and B, where we discuss the line iden- 
tification and give the dipole moment for the different 
observed vibrational transitions. The continuum emis- 
sion is analyzed in §4. The results concerning C2H2 and 
its isotopologues are presented in §5 and §6, and those 
related to HCN and H 13 CN in §7. The results obtained 
for C2H2 and HCN, their uncertainties, and the sensitiv- 
ity of the model to the physical conditions of the CSE 
are discussed and analyzed in §8. Finally, the results are 
summarized in §9. 

2. OBSERVATIONS 

We observed IRC+10216 with TEXES, the Texas 
Echelon-cross-Echelle Spectrograph (Lacy et al. 2002), 
at the NASA Infrared Telescope Facility on 2002 Dec 12 
(UT), corresponding to an IR stellar phase </>ir, ~ 0.08 
(following Monnier et al. 1998). We used the TEXES 
high resolution echelon grating with a first order grating 
as the cross-disperser. In this mode, we obtained a spec- 
tral coverage of roughly 0.25 /im per setting. To cover 
the entire range shown here, 11.6 to 13.9 /xm, we required 
10 separate settings and stored the data in 10 different 
files. 

At these wavelengths, the echelon orders are larger 
than our detector, meaning that there are gaps in the 
spectrum between orders. The gaps increase in size 
toward longer wavelengths. Telluric features become 
stronger and more frequent toward 13.9 fim. In regions 
where the telluric atmosphere was nearly opaque (<5% 
transmission) we discarded the data, resulting in addi- 
tional gaps. 

When using the first order grating as the cross dis- 
perse^ the TEXES slit must be very small to prevent 
orders from overlapping. For these observations the slit 
was 2'.'2 long. With the short slit, we nod IRC+10216 
off the slit for sky subtraction. Sky conditions were good 



enough that we do not believe this introduced any sig- 
nificant systematics. 

The data were reduced with the standard TEXES 
pipeline (Lacy et al. 2002). The pipeline corrects opti- 
cal distortions, combines nod pairs, removes spikes, flat- 
fields, performs correction for telluric absorption, estab- 
lishes a frequency scale using telluric features, and ex- 
tracts a spectrum. We normalized each spectrum before 
analysis using a fourth-degree polynomial to estimate the 
baseline. 

Typically, observations of asteroids provide the best 
telluric corrections. However, IRC+10216 is so much 
brighter than the available asteroids that we would have 
limited our signal-to-noise ratio by using asteroid mea- 
surements. Instead, we used a blackbody-sky difference 
spectrum to correct for the atmosphere. In the ideal case 
of the atmosphere, telescope, and blackbody being at a 
uniform temperature, the difference spectrum indicates 
the absorption from the sky. In our experience, this pro- 
cedure for atmospheric correction works fairly well except 
for telluric water vapor, which can change on a short time 
scale and from one line of sight to another. 

The observed spectrum covers the wavenumber range 
720 to 864 cm -1 . Taking into account the molecules 
detected in this source (Keady & Ridgway 1993; Cer- 
nicharo et al. 2000), there are several candidates that 
could contribute to the observed features: C2H2, HCN, 
and SiS, including their isotopologues. All of them have 
vibrational bands in the observed wavelength range (see 
Figures 1—6). Although many lines of SiS and other 
molecules have been identified in the Figures, they will 
be analyzed and studied in a forthcoming paper. 

The spectra was corrected from the source movement 
by identifying and modeling most of the C2H2 i^Re and 
HCN i^Re lines in the spectrum, and calculating the 
mean shift between the observed frequencies and those 
published from laboratory work (Rothman et al. 2003). 
This shift is ~ 45 km s" 1 . However, because of small 
uncertainties in the observational process, the data in 
a few files needed for an extra correction smaller than 
2 km s _1 . Moreover, we had to correct the central fre- 
quency of some C2H2 and HCN bands taken from the 
HITRAN Database by blue-shifting them in less than 
1.5 km s _1 . A sample of the velocity-corrected data is 
shown in Table 1. 

In order to fit the continuum, we have used ISO/SWS 
observations carried out on 1996 May 31 (UT), that cor- 
responds to 0ir ~ 0.34. 

3. THE MODEL 

3.1. The Physical Structure of IRC+10216 from 
Previous Work 

IRC+10216 is ejecting matter (gas and dust) at a rate 
of 1 - 2 x 10~ 5 Mq yr" 1 (Keady et al. 1988; Cernicharo 
et al. 1996b, 1999). The derived dust ejection rate is 
between ~ 2 - 4 x 10~ 7 M yr _1 (Ridgway & Keady 
1988; Men'shchikov et al. 2001). The star is pulsating 
with a period of 636 ± 3 days (Ridgway & Keady 1988; 
Dyck et al. 1991; Jones et al. 1989). The stellar effective 
temperature, T c g, quoted in the literature varies from 
author to author: 2330+350 K (Ridgway & Keady 1988), 
2200 ± 150 K (Ivezic & Elitzur 1996), 1915 K at phase 
0.16 and 2105 K at phase 0.27 (Bergeat et al. 2001), and 
2800 and 2500 K for maximum and minimum brightness 
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Fig. 1. — The observed mid-infrared spectrum of IRC+10216 in the range 720-744 cm -1 with resolution ~ 10 — 2 cm -1 (~ 3 — 4 km 
s" 1 ). The most important transitions of C 2 H 2 (red), HCN (blue), H 13 CCH (green), and H 13 CN (light blue) are indicated. The sensitivity 
is high enough to detect molecular species with abundances as low as 10~ 8 — 10~ 7 relative to H2. Gaps occur between grating orders and 
in regions of high atmospheric opacity. Many of the unlabeled lines remain unidentified. 
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Fig. 4. — The observed mid-infrared spectrum of IRC+10216 in the range 792-816 cm with resolution ~ 10 cm" _1 (~ 3-4 km s" 1 ). 
The most important transitions of C 2 H 2 (red), HCN (blue), H 13 CCH (green), and H 13 CN (light blue) are indicated. The -y-axis scale has 
been changed beyond 810 cm -1 to show the molecular features more clearly. 
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Fig. 6. — The observed mid-infrared spectrum of IRC+10216 in the range 840-864 cm~' with resolution ~ 10 cm" - 1 (~ 3-4 km s" 1 ). 
The most important transitions of C2H2 (red), HCN (blue), H 13 CCH (green), and H 13 CN (light blue) are indicated. The -y-axis scale has 
been changed with respect to Figures 1-3 to show the molecular features more clearly. 
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TABLE 1 

Observed Spectrum towards IRC+10216 



1 N 

Freq. (cm ) 


Normalized Flux 


— ; 

Atmospheric Trans. 


Model 


721.579 


1.023 


0.1810 


1.000 


721.581 


1.089 


0.1877 


1.000 


721.583 


1.126 


0.1940 


1.000 


721.586 


1.121 


0.2009 


1.000 


721.588 


1.111 


0.2068 


1.000 


721.590 


1.079 


0.2127 


1.000 


721.592 


1.019 


0.2166 


1.000 


721.594 


0.9778 


0.2159 


1.000 


721.597 


0.8599 


0.2144 


1.000 


721.599 


0.7520 


0.2155 


1.000 



Note. — Observed spectrum towards IRC+10216 ranging from 
721 to 864 cm -1 . Column 1 contains the observed frequencies in 
cm -1 corrected from velocity shifts due to proper motions of the 
source, frequency calibration and uncertainties related to the band 
center of laboratory determined frequencies. The observed flux, hav- 
ing removed the baseline, can be found in Column 2. Column 3 
accounts for an estimation of the atmospheric transmission during 
the observations at the considered frequencies. The model results 
have been included in Column 4 to allow comparisons. See the text 
for details. [The complete version of this table is in the electronic 
edition of the Journal. The printed edition contains only a sample.] 

respectively (Men'shchikov et al. 2001). 

The distance of IRC+10216 is poorly established; val- 
ues vary between 120 and 300 pc (Doty & Leung 1997; 
Keady et al. 1988; Weigelt et al. 2002; Bergeat et al. 2001; 
Loup et al. 1993; Cernicharo et al. 2000; Hcrbig & Zap- 
pala 1970). From intcrfcrometric and lunar occultation 
data, Ridgway & Keady (1988) have obtained a stellar 
angular radius of 0'.'019 ± 0'.'003, Keady et al. (1988) as- 
sumed 0'.'023, the value derived by Monnier et al. (2000a) 
is 0'.'022, and Men'shchikov et al. (2001) suggest an an- 
gular stellar radius ranging from 0'.'014 to 0'.'018 over the 
whole period of pulsation (corresponding to ~ 815 R for 
a distance of 200 pc derived from observations at 1.65, 
2.2, 3.15, and 4.95 /im; ~ 970 R Q with a distance of 
200 pc; ~ 635 R Q for a distance of 135 pc from 8-12 /im 
observations; and a radius ranging from 390 to 500 R Q 
with an assumed distance of 130 pc, from measurements 
between 0.6 ^m and 6 mm, respectively). 

Infrared observations during a lunar occultation (Ridg- 
way & Keady 1988) showed that the dusty inner envelope 
of IRC+10216 is asymmetric. There is a large amount of 
dust in the equatorial plane and two bright lobes along 
the poles where the dust density and resulting extinction 
are lower. Interferometric observations by Weigelt et al. 
(1998) revealed at least four clumps in the lobes. Later 
observations showed that these clumps evolve with time 
(Weigelt et al. 2002; Tuthill et al. 2005) with timescales 
of ~ 1 year (Men'shchikov et al. 2002). The large scale 
density profile of the dust reveals nearly concentric shells 
corresponding to increased ejection of matter over peri- 
ods of 200 - 800 years lasting 20 - 40 years each (Mau- 
ron & Huggins 1999, 2000; Murakawa et al. 2002). The 
latter ejection episodes are compatible with results by 
Men'shchikov et al. (2001) that suggest that the star has 
experienced at least two episodes of high mass loss over 
the last 1000 years. 

From the work of Keady et al. (1988), dust grains are 
formed by amorphous carbon (AC) with some inclusions 
of SiC and perhaps other components containing Mg or 
S. SiC condenses in the photosphere, forming the seeds 
of dust grains (the condensation temperature of SiC is 



~ 2000 K - Men'shchikov et al. 2001). Ivczic & Elitzur 
(1996) have calculated the relative amounts of the main 
components of the dust grains to be ~ 95% amorphous 
carbon, 3 - 8% SiC, and less than 10% MgS. Other au- 
thors have suggested that the molecular bands detected 
in the continuum can be produced by molecules more 
complex than MgS (e.g., Men'shchikov et al. 2001). 

The gas is accelerated to ~ 2 km s _1 near the photo- 
sphere (innermost envelope). At some distance from the 
star the carbonate material condenses and the gas is ac- 
celerated again until it reaches a velocity of ~ 11 km s _1 . 
Ridgway & Keady (1988) have obtained an inner ra- 
dius for the condensation of carbonate material of 5 R* . 
At ~ 11 R* refractory Mg and/or S-bearing molecular 
species condense and the gas reaches an expansion ve- 
locity of ~ 14 km s _1 . Keady et al. (1988) propose a 
turbulence in the innermost CSE of ~ 5 km s _1 and a 
terminal turbulence velocity of ~ 1.0 km s _1 . Lower val- 
ues for the terminal turbulence velocity can be found in 
the literature, e.g., 0.9 km s _1 (Huggins & Healy 1986) 
and 0.65 km s" 1 (Skinner et al. 1999). 

3.2. Model Description 

A good approach to the spectral line radiative transfer 
problem in a CSE is given by the LVG method. It is 
very fast and produces good results but assumes that the 
linewidth of a line is negligible compared to the velocity 
gradient between very close emitting regions of the CSE. 
Hence, they cannot be used to accurately model warm 
regions like the innermost CSE (T K ~ 1000 - 2500 K) 
where the linewidth can be larger than the expansion 
velocity gradient. 

The most exact model we can develop should be a non- 
local model which involves both the numerical resolu- 
tion of the statistical equilibrium and radiative transfer 
equations. Studies at radio wavelengths have success- 
fully applied this method to CSEs (ALI codes by, e.g., 
Justtanont et al. 2005 and Monte Carlo codes by, for ex- 
ample, Gonzalez-Alfonso & Cernicharo 1997; Crosas & 
Menten 1997; Schoicr & Olofsson 2006), although they 
have needed a significant amount of CPU time even when 
only a few vibrational and rotational levels have been 
considered. Nevertheless, the most important reason 
that led us to reject the use of a non-LTE code is the lack 
of ro- vibrational collisional rates quoted in the literature 
for C2H2 and HCN. We have identified ro-vibrational 
lines of C2H2 and HCN involving vibrational levels with 
energies up to 2200 cm -1 that are created close to the 
star. In the innermost CSE, the kinetic temperature is 
high enough to significantly populate these vibrational 
levels having energies between 3000 and 4000 cm -1 or 
more. For example, C2H2 has about 50 vibrational lev- 
els below 4000 cm -1 . Therefore, a thorough model of 
the envelope should account for many vibrational and 
rotational states. The large number of ro-vibrational lev- 
els to be considered, makes the problem computationally 
impractical with these kind of methods. We can signifi- 
cantly reduce the computing time by avoiding exact so- 
lutions of the statistical equilibrium equations. We do 
this by assuming a temperature dependence with radius 
(supported by previous work or physical considerations 
and by allowing several parameters free in order to fit 
the observed lines. A reasonably good estimation to the 
excitation temperatures could be achieved taking advan- 
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TABLE 2 

Identified and Modeled Lines of C2H2, H 13 CCH, HCN, and H 13 CN 



Ro- Vibrational Band 
Upper level Lower level Branch 



# Detected Lines 
C2H2 H 13 CCH 



Ro- Vibrational Band 
Upper level Lower level Branch 



# Detected Lines 
C2H2 H 13 CCH 



^5 (^u) 


lj.O.(CTg ) 




OU I^OU ) 


oy 


.ZD ) 






H f 


ACl ( C )A\ 






Ve 


91 ( (W 


1 7 
1 1 


' (Y\ 
\ U J 






p* 

^7 


1 1 C\\ 
1 { UJ 








fe 


2 ( 0) 


1 


v °) 




2^4 (ct J) 


Re 


7 ( 0) 


( 0) 




Va(t{ n) 


R e 


46 (15) 


11 


' 9) 






Re 


4 ( 0) 


( 0) 


U4 + V5 {&u ) 


U 4 {lTg) 


Rf 


42 (25) 


9 


: 4) 






«/ 


2 ( 0) 


( 0) 






Qe 


20 ( 0) 


2 


' 0) 


2^/1 -1- L'c ft^i, 1 


2va(8„) 


Re 


7 f 0) 


( 0) 






Pf 
J 


3 ( 0) 





' 0) 






Qe 


2 ( 0) 


( 0) 




Va(tV ) 


Re 


37 (25) 


9 


' 7) 






Rf 

J 


11 ( 0) 


( 0) 






Qf 


3( 0) 





; 0) 




va + v 5 (at) 


Qf 


2 ( 0) 


( 0) 






Pe 


2( 0) 


1 


; 0) 


i/4 + 2^-5 (7T 9 ) 


Re 


2 ( 0) 


( 0) 






Rf 


44 (27) 


6 


: 5) 


i/4 + 2v 5 (TT g ) 


VA + 1/5(0-";) 


Rf 


5 ( 0) 


( 0) 






Qe 


13 ( 0) 


1 


: 0) 


I/ 4 + 2f 5 (0g) 


1/4 + 1/5 (<5u) 


Re 


7( 0) 


( 0) 


2v b (<t+) 




Pj 


1 ( 0) 


1 


; 0) 






Rf 


8 ( 0) 


( 0) 


Vh(~Kv) 


Re 


34 (18) 


2 


: 2) 


3j/5(tt u ) 


2i/ 8 (^) 


Re 


2 ( 0) 


( 0) 






Q f 


3( 0) 





: 0) 


3i/5(0u) 


2v 5 {5 9 ) 


Re 


5 ( 0) 


( 0) 


2v 5 (8 a ) 




Re 


31 (16) 


5 


: 5) 






Rf 


4( 0) 


( 0) 






Qe 


2( 0) 





: 0) 












Upper level 


Lower level 


Branch 


HCN 


H 13 


CN 


Upper level 


Lower level 


Branch 


HCN 


H 13 CN 


1/2(77) 


G.S.(<7+) 


Re 


37 (19) 


7 


: 5) 


2lAj((J) 


V2 (t) 


Re 


23 (14) 


( 0) 






Qe 


3( 0) 





; 0) 






Rf 


22 (17) 


( 0) 


2l/ 2 ((T+) 


1/2(1") 


Re 


11 (10) 


( 0) 






Qe 


1 ( 0) 


( 0) 



Note. — The numbers without parentheses are the detected lines while the numbers within parentheses are the modeled 
lines. Many lines for each band are contaminated by blending with other features, overlapping with telluric features, and 
other effects such as echelon order boundaries. The detected lines involving high-J and high energy ro-vibrational levels will 
be published in a forthcoming paper. Those chosen to be fitted are the best quality lines for each band. The spectroscopic 
notation used is described in Appendix A. 



tage of the large number of observed ro-vibrational lines, 
allowing us to estimate also the populations of many ro- 
vibrational levels below 2500 cm -1 . This handmade pro- 
cessing technique accounts for all the physical phenom- 
ena having an influence in the molecular level popula- 
tions at the nearby environment of each position in the 
envelope, for example, the collisional rates. 

Concerning the geometry of the envelope, modeling 
the important deviations from spherical symmetry in the 
innermost CSE requires a large number of parameters 
which have limited impact in our fits due to the low angu- 
lar resolution of our observations. Nonetheless, the likely 
complex velocity field of the gas might introduce some 
small features in a well defined region of the emission 
component of the the high-excitation line profiles. Un- 
fortunately, the signal-to-noise ratio is not large enough 
to unmistakcably identify the signs of the inner structure. 
Moreover, the envelope approaches spherical geometry at 
large scales. Therefore we have ignored the complex in- 
ner dust structure and assumed spherical symmetry for 
the whole CSE. 

These approximations allow us to solve the problem 
with one dimensional calculations. We choose a two di- 
mensional coordinate system with the y axis parallel to 
the line of sight pointing away from the Earth and the x 
axis perpendicular to y in an arbitrary direction. Then, 
we study the evolution of just one ray of light parallel to 
y for each x and build the final spectrum using the sym- 
metry of the envelope. The emerging intensity at each 
x is calculated by summing the emission of each region 
backwards from the star, so that the radiation emitted 
by a region is affected by the optical depth of those in the 
foreground. Due to its importance in the mid-IR, dust 



is included in the calculations and radiatively coupled to 
the gas. The last step is to multiply the emerging inten- 
sity by the point spread function (PSF) of the telescope, 
take into account the slit dimensions of the spectrometer 
and to convolve the resultant flux with the frequency re- 
sponse of the detector. By virtue of simplicity, we have 
adopted a gaussian profile for both the PSF and the de- 
tector response with the HPBW of the telescope and the 
width of each channel of the detector as FWHM, respec- 
tively. 

We have adopted the velocity field structure proposed 
by Keady et al. (1988) which consists of three different 
zones: Region I starts at the stellar photosphere and 
ends at the distance of the first dust formation layer, 
Rdi; Region II starts at that radius and extends to the 
position of the second dust formation layer, Rd2; finally, 
the rest of the envelope is Region III. We have initially 
assumed that the two dust formation regions are located 
at 5 and 15 R* and have a thickness of 1 R* (Keady et 
al. 1988). We allow these values to vary in the model 
to get the best fit for the dust emission (continuum of 
the star) and the molecular features. We have divided 
the envelope into a large number of concentric shells to 
follow the rapid variation of physical conditions. Using 
a logarithmic step for the radius increment, most shells 
are placed in the inner and middle regions (Regions I 
and II), where the temperature and density gradients 
are more important. 

We have tried to find some simple laws for the vari- 
ation of molecular abundances, density, and gas kinetic 
and dust grain black-body temperatures as a function 
of r that best reproduce the observed continuum and 
line intensities and profiles. We have assumed that all 
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the temperatures (dust, kinetic, vibrational, and rota- 
tional) follow a continuous radial dependence r~ a , where 
a could be different for each region of the envelope and 
for each ro- vibrational level (a — > a v j), depending on 
the considered temperature. In the case of the dust 
temperature, we have assumed that a remains constant 
over the whole dusty CSE. The parameter a (for Re- 
gions I and II, and for each temperature except the dust 
temperature) is completely determined through the in- 
put of the corresponding temperature at R» and Rdi 
for Region I, and Rdi and Rd2 for Region II. In Region 
I T x (R*) /T x (Rdi) = (R^/Rdi)-"- 1 and in Region II 

T» (Rdi) /T x (Rda) = (Rdi/Rd 2 )" Q - n , where x can be 
dust, kinetic, vibrational, or rotational. However, a in 
Region III cannot be obtained accurately from a model- 
ing of different lines. In consequence, we have assumed 
a = 1.0 for all the temperatures in Region III (Doty & 
Leung 1997). As we will see below, this hypothesis is 
compatible with the fits. On the other hand, the gas 
density is assumed to satisfy the continuity equation: 

M = 47rr 2 w oxp ^TOH 2 ^H 2 (r) (1) 

where p, = ^] miXi/mn 2 , mi is the mass of the i th most 
abundant species and x\ its abundance with respect to 
H2. Since the most abundant species after H2 are He 
and CO (with x ~ 0.2 - solar abundance, Cox 2000 - 
and ~8x 10~ 4 , respectively) then, ~ 1.4. 

For dust grains, we have assumed a static density pro- 
file following a r~ 2 variation law: 

Mr) = ^(^-) 2 (2) 
a\ Rdi \ r J 

where t\ and a\ are the optical depth and the absorption 
of a dust grain at wavelength A respectively, Ao is a fixed 
wavelength and t\ is an input parameter of the model. 
The wavelength Ao that we have adopted is 11 /im. The 
dust opacity is derived from the optical properties of 
amorphous carbon (AC) and silicon carbide (SiC) at each 
wavelength and the emission is computed from the dust 
opacity and temperature in each volume clement. 

The emission and absorption in an elementary integra- 
tion step is calculated from the adopted velocity field, the 
corresponding vibrational and rotational temperatures, 
the H2 density and the molecular abundances. For the 
gas, the adopted line frequencies used for the line iden- 
tification are given in Appendix A, while all the data 
relevant to line intensities (opacity, dipole moment, and 
partition function) are in Appendix B. 

The line profile is assumed to be gaussian with a ther- 
mal linewidth over the external CSE because the dust 
and the gas are expected to expand at their terminal 
velocity in a steady flow. In contrast, the linewidth 
next to the photosphere has not been well determined 
by previous observations of IRC+10216. The matter 
ejection near the star could increase the microturbu- 
lent velocity adding a non-thermal contribution to the 
linewidth. In our code, this non-thermal contribution 
has been implemented as At! m ,ie _ ' r ~ 1 )'^, where Az; mj i 
is the non-thermal microturbulent velocity over the pho- 
tosphere and £ is a characteristic length (r and I are 
measured in R*). Since the lines affected by tempera- 
ture and/or microturbulence present gaussian profiles, 
the resultant linewidth is (Aw) 2 = (Ai>temperaturc) 2 + 
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Fig. 7. — Optical depth per unit length integrated over frequency 
for several lines of the i/5(7r u )R e branch of C2H2- The optical depth 
for low- J lines is still large at 50 R* . However, high- J lines (J > 20) 
will be mainly formed in the inner and middle regions (Region I, 
72-j, and Region II, 72jj)- Hence, we can obtain information about 
the outer envelope (Region III, Tim) from low-J lines while high- 
J lines are sensitive to the physical and chemical conditions near 
the star. Choosing the correct line can help improve the accuracy 
of the determination of the parameters in a particular region (see 
Table 7). The thick red line represents the velocity field. The sharp 
increase in opacity at ~ 5 R* corresponds to the change in C2H2 
abundance in Region I found in this work. 

(Aw microturbuloncc ) 2 . In order to obtain physically signif- 
icant parameters from the code, the parameter we have 
used is the total linewidth in the photosphere, Av±. High 
energy ro- vibrational lines, which are formed at the pho- 
tosphere, could carry some information on these param- 
eters. 

As many observed lines are modified by many weak fea- 
tures, telluric effects, and overlaps with other lines, the 
fits derived by the minimization of x 2 would be worse 
than those obtained here. In addition, the model de- 
pends on many parameters and minimizing the x 2 func- 
tion for all the selected lines is computationally unrealis- 
tic. In fact, we were able to obtain reasonably good eye 
fits which are in agreement with those derived by mini- 
mizing the \ 2 function when diminishing the number of 
parameters (as a first approximation). 

However, in the case of the uncertainties of the param- 
eters (see §8), we have used a numerical method. Given 
an observed line, several parameters could similarly affect 
the synthetic profile masking their effects and hindering 
an accurate determination of all their values that reason- 
ably fit the line. Thus, we could calculate the maximum 
value adopted by a given parameter, by considering it as 
a function of the others (having imposed a restriction on 
the synthetic profile with respect to the observed one; for 
example, setting the function x 2 to be equal to a given 
number or forcing the synthetic profile to differ from the 
observed one in less than a given quantity), and by fol- 
lowing the gradient of this function. It is also possible 
to find the minimum value of that parameter by follow- 
ing the field opposite to the gradient. This is a good 
method to be used for a large number of parameters and 
for a few lines, since it is possible to approach the maxi- 
mum (minimum) in several iterations avoiding unneces- 
sary calculations and spending a reasonable amount of 
CPU time. 

In our models we have adopted a stellar temperature 
of 2330 K, a distance of 180 pc, an angular radius for 
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the star of 0'.'019 (corresponding to a stellar radius of 
735 Rq) and a mass loss rate of 2.1 x 10~ 5 M© yr _1 . 

The region between 1 and 100 R* is shielded against 
Galactic UV photons by the outer dust envelope. Hence, 
we have considered the radiation field coming only from 
the star and the dust component of the CSE. 

Although the envelope can extend to more than 
1000 R», the dust in the outermost part of the CSE 
does not contribute significantly to the continuum flux 
in the mid-infrared (an envelope with a maximum radius 
of 600 R* divided in 300 shells is enough to find a good 
fit for the continuum). 

The vibrational and rotational temperatures decrease 
so fast that only the emission and absorption of the low- J 
levels of the ground vibrational state of C2H2 and HCN 
could be affected by the choice of the external radius in 
our models (e.g., at 100 R* the kinetic temperature is 
~ 85 K and the most populated level of C2H2 is J ~ 
4 — 5, sec Figure 7). The outer radius of the CSE for 
fitting the lines has been fixed to 300 R* since there is no 
significant modification due to the CSE region between 
300 and 600 R». A total of 100 layers were modeled. 

The absorbing feature at terminal velocities is mainly 
produced in Region III for low- J ro- vibrational lines. Op- 
tically thin absorbing lines are produced essentially in 
Region II, related to the (—10, 10) km s _1 velocity range 
in the profiles, setting the absorption maxima at veloci- 
ties around —8 and —10 km s _1 . The innermost region of 
the CSE (Region I), and the fact that the dust continuum 
is largely formed outside of this radius, has little effect in 
low- J line profiles due to its small angular size compared 
with the rest of the CSE. The kinetic temperature and 
density are high enough to assure LTE throughout the 
envelope for J < 20 within vibrational levels up to 2v§ 
for C2H2 and up to 2j/ 2 for HCN. The validity of this 
assumption has been checked in our models and found 
to be acceptable. However, high energy ro-vibrational 
levels are populated almost completely in Region I, al- 
lowing us to derive physical parameters for the regions 
closest to the star. Consequently, the sensitivity of the 
observed line profiles to the physical conditions of the 
gas and the large variation of these parameters across 
the envelope, limit the radial resolution we can achieve 
for the molecular abundance profiles. 

4. CONTINUUM EMISSION 

The dust properties affect considerably the molecular 
excitation due to radiative coupling between the dust and 
the gas. P-Cygni profiles arise in all the lines created in 
the inner CSE. In particular, the C2H2 and HCN lines 
show this kind of lineshape (see Figures 1-6). The line 
profile depends on the dust parameters such as the dust 
temperature, composition, size of the dust grains, and 
absorption and scattering cross-sections. 

For the wavelength range under consideration (MIR, 
A ~ 11 — 14 (um), the size of the grains has a small ef- 
fect on the continuum. This parameter can be neglected 
at larger wavelengths (FIR and radio) but is much more 
important in the NIR range (ivezic & Elitzur 1996), i.e., 
A < 1 — 2 fim. We have assumed that dust grains are 
spheres with a constant radius not larger than 0.1 /zm and 
calculated their opacity using Mie Theory (e.g., Hoyle 
& Wickramasinghc 1991). According to the theory, the 
scattering cross-section is small compared to the absorp- 



tion cross-section for dust grains with diameter small 
compared to the wavelength of radiation. Since we are 
considering A ~ 11 — 14 //m radiation, we ignore scatter- 
ing. The composition of the dust grains are guessed to 
be amorphous carbon (AC) and silicon carbide (SiC) (see 
§3.1). The complex refractive index of AC as a function 
of A has been taken from Rouleau & Martin (1991). Lab- 
oratory works on SiC (Mutschke et al. 1999) suggest that 
SiC is probably crystalline in space. Unfortunately, the 
difficulties in measuring the optical constants for crys- 
talline SiC have led us to adopt the laboratory data for 
the refractive index of amorphous SiC in our models. For 
the dust temperature, Td, we initially assume the r~ 0A 
dependence, proposed by Ridgway & Keady (1988). Ac- 
tually, the Td dependence on r is steeper in the innermost 
CSE than at greater distances because the opacity of the 
dust grains is larger at higher frequencies, i.e., closer to 
the star. We have adopted the power law shown above 
as an approximation. The exponent is a very sensitive 
parameter and can be determined with a high accuracy 
(see §8). Finally, concerning the dust optical depth, t\, 
the latter authors derived a value rx — 1 at 11 /im, Ivezic 
& Elitzur (1996) obtained r A ~ 0.32-0.40 at 10 /im, and 
Monnier et al. (2000a), obtained t\ ~ 0.66 at 11.15 /urn. 

In this work, we have fitted the continuum data of 
IRC+10216 obtained by ISO/SWS, adjusting the dust 
parameters as described above. The ISO observations 
used for that purpose are those quoted by Cernicharo et 
al. (1996b, 1999). We have assumed that the difference 
between the stellar phases of the observations (see §2) 
with ISO/SWS, and IRTF/TEXES, does not introduce 
any important effect on dust and molecular properties, 
and that the derived Td from the ISO data is similar to 
those prevailing at the moment of the TEXES observa- 
tions. 

We have fitted the continuum between 7 and 27 /zm. 
Fitting to shorter wavelength ISO data is inappropriate 
because scattering will be important and the model does 
not include scattering. The brightness of the central star 
has a small impact on the observed continuum. The max- 
imum contribution to the continuum from the star is at 
5.31 /im (F v ~ 1075 Jy) and the emission at 10 /im is 
~ 1.73 times lower. A small modification of the stellar 
temperature does not significantly affect the flux emerg- 
ing from the source. Note, however, that a large change 
in the stellar temperature could modify the physical and 
chemical properties of the envelope and, hence, of the 
emergent flux. In any case, it seems obvious that the 
observed mid- and far-infrared emission in IRC+10216 
comes mainly from the dusty envelope. 

We obtain from our analysis that the inner and outer 
dust formation shells are located at radii 5.2+q g R* and 
21 ± 3 R* (O'.T and 0'.'4), respectively (see Figure 7, Ta- 
ble 7, and §8). A value, Rdi = 5.2 R* is derived from 
fitting only the continuum. However, Rd2 can be deter- 
mined with more or less precision through fitting several 
molecular lines, while the continuum of the star displays 
little information about this outer dust formation shell. 
Studying how changes to the dust density, n d , affect Rd2 
led us to very inaccurate results because the variation of 
rid over this shell only slightly affects the continuum. 

Several studies have determined Rdi since the discov- 
ery of IRC+10216. Keady et al. (1988) proposed that 
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Rdi ~ 3 R» equivalent to ~ 0'.'057, with a stellar ra- 
dius of ~ 970 R© and a distance to the star of 200 pc. 
Later, Ridgway & Keady (1988) found Rdi ~ 5 R* and 
an angular stellar radius of ~ 0'.'019, implying an angular 
radius of the inner dust formation zone, ctyi, of ~ 0'.'095. 
Monnier et al. (2000a) found ayi — 0'.'15 at a distance of 
135 pc and a* ~ 0'.'022, meaning that Rdi ~ 6.8 R©, and 
Ivezic & Elitzur (1996) determined ctyi to be ~ 0'.'22 at 
maximum luminosity and ~ 0'.'15 at minimum luminos- 
ity. Consequently, our result is in good agreement with 
the values quoted in the literature. For Rd2, the angu- 
lar radius of the outer dust formation zone, Keady et al. 
(1988) found R d2 ^ 14 R« (a d2 ~ 0'.'32). We derived 
Q^d2 — 0'.'40, a larger value than found by Keady et al. 
(1988) and Keady & Ridgway (1993). 

The optical depth derived with our model at 11 /im is 
0.7, similar to that proposed by Monnier et al. (2000a), 
t(A = 11.15 /im) = 0.66, and larger than the value pro- 
posed by Ivezic & Elitzur (1996), who derived t(A = 
10 //m) ~ 0.32 — 0.40. On the other hand, the value for 
the temperature of the innermost dust formation shell, 
T d i, derived in this work is 850±25 K (see §8 and Table 7 
for a discussion on the errors), lower than that obtained 
by Ridgway & Keady (1988) (T dl = 1040 ± 100 K) and 
Groenewegen (1997) (T dl = 1075±50 K), similar to Mon- 
nier et al. (2000a) (T d i = 860 K), and higher than Ivezic 
& Elitzur (1996) (T A1 = 750 ± 50 K). The exponent of 
the Td law is derived to be 0.39, quite similar to that 
proposed by Keady et al. (1988). 

The best fit to the dust composition is 95% amorphous 
carbon and 5% amorphous silicon carbide. Modifying the 
proportion of SiC changes the ratio of the predicted flux 
at 11 /im, with respect to the rest of the continuum. Nev- 
ertheless, this change in composition could be balanced 
by modifying the optical depth and dust temperature at 
Rdi. Fitting the whole continuum is necessary to get 
more realistic percentages. 

Taking into account the different observing periods and 
the wavelength coverage in these different data sets, we 
consider that our estimate of Td and opacity in Rdi are 
representative of the physical conditions in IRC+10216. 

5. C 2 H 2 , MODELING 

In order to obtain the abundance of C2H2 and the 
physical parameters of the CSE, we have selected several 
sets of C2H2 lines that could be sensitive to the derived 
parameters (see Table 2 and Figure 7) . The assumed gas 
velocity profiles and computed line profiles arc in very 
good agreement (see Figure 8). The C2H2 abundance 
and physical conditions from the fits are listed in the 
following subsections. 

5.1. Velocity Profile 

The simultaneous modeling of different J ro- 
vibrational transitions of the same band requires a three- 
region velocity field (Keady et al. 1988). As can be seen 
in Figure 7, lines with higher energy are formed closer 
to the photosphere. Consequently, the effect of the ex- 
pansion on the lineshapc changes with J. In many cases, 
a variation of the velocity in one of these regions mod- 
ifies the lineshapes significantly. For example, the gas 
velocities in Regions II and III control the position of 
the intensity minimum, so that lines formed closer to 
the star represent the bulk absorption at lower velocities 



than the lines formed at larger distances. 

On the other hand, the model is very insensitive to 
variations in the thickness of the acceleration shells (see 
§8). The linewidths in these regions are similar to the 
velocity-increasing in the acceleration zones, and hence 
hide the acceleration effects on the profiles. The small 
spatial extent of the innermost region and the limited 
angular resolution of our data, prohibit a better deter- 
mination of the physical and chemical conditions in Re- 
gion I. We have found that the best fit to the lines can be 
achieved with the following velocity profile (see Figure 7): 

{5 1 < r/R« < 4.7 

5-11* 4.7<r/R»<5.7 

11 5.7 < r/R* < 20.75 (3) 

11 - 14.5* 20.75 < r/R* < 21.75 

14.5 21.75 < r/R* 

where * means constant velocity gradient, dv e /dr between 
the inner and outer radii of the corresponding zone. 

5.2. Linewidth 

When microturbulence is included in the calculations, 
the emission from Region I spreads and the intensity of 
the lines reaches its maximum at higher velocities (to- 
wards lower frequencies), compared with the fits with- 
out microturbulence. Unfortunately, fits for low- and 
intermediate- J ro- vibrational levels are relatively insen- 
sitive to the total linewidth in the photosphere, Avi, 
and the characteristic length, £. A substantial change 
in the fitted profiles requires unrealistic values for these 
parameters. Reasonable values, for example, are Avi = 

5 km s _1 and I = 1.5 R*. The value of I is acceptable 
due to the physical phenomena thought to occur near the 
photosphere, i.e., a moderate pulsation of the central star 
and the appearance of sound waves (Bowen 1988; Pijpers 

6 Hearn 1989a; Pijpers & Habing 1989b). However, the 
fits obtained with these parameters are not more accu- 
rate than those calculated without microturbulence. In 
fact, when choosing A«i = 30 km s _1 and I = 1.5 R„, the 
emission can be fitted more accurately. We will discuss 
this further in §8. 

5.3. Kinetic Temperature, Tk 

C2H2 low and intermediate energy levels are highly 
populated in each vibrational level, and many of the ro- 
vibrational transitions in which they are involved are op- 
tically thick and therefore hide information (e.g., low- 
J transitions in the fundamentals in Regions I and II). 
However, there are still many lines in LTE that are not 
optically thick in the hot bands over the whole enve- 
lope and in the fundamental band in Regions II and 
III (see Figure 7). We can determine the kinetic tem- 
perature profile by fitting simultaneously several low- 
and intermediate- J lines from the fundamental band and 
relatively low energy hot bands. The relationship be- 
tween the intensities of adjacent ro- vibrational lines in 
LTE within a given band allows us to derive Tk immedi- 
ately. In doing so, it is necessary to derive the vibrational 
temperatures and the molecular abundances at the same 
time as Tk. We will discuss this determination in §5.5. 

We have fixed the kinetic temperature at r = 1 R* 
equal to the temperature of the central star. For Regions 
I and II, Tk at Rdi and Rd2 are free parameters that 
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Fig. 8. — Fits of some lines of C2H2. The black lines are the observed spectra while the red lines are the fits. The data have been 
corrected from observational frequency deviations adding/suppressing the mean difference between the observational and the laboratory 
frequency of many lines (see §2). The existing discrepancies are small and can be assumed as blends with other lines and/or observational 
errors (the spectroscopic notation is described in Appendix A). 
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completely determine the value of the exponent. The 
temperature exponent obtained is 0.58 for both Regions 
I and II, while the fits support the assumption of a = 1.0 
over Region III. 

The kinetic temperature profile derived from the best 
fit to all C2H2 lines is: 

( 2330 (r/R»)-°- 58 Region I 
T K (r) = \ 900 (r/R dl )-°- 58 Region II (4) 
[ 400 (r/R d2 )- 100 Region III 

5.4. Rotational Temperatures, T rot 

An accurate estimate of the rotational temperatures 
requires a sophisticated model including molecular pa- 
rameters, such as collisional rates between ro- vibrational 
levels, that are currently unavailable. Therefore, we have 
used our data to derive a first estimate of rotational tem- 
peratures for all ro-vibrational transitions in the three 
zones defined above. The procedure to derive T rot pro- 
files is as follows: 

1. Tk and the vibrational temperatures are derived by 
modeling the ro-vibrational transitions that involve 
both low and intermediate energy rotational levels 
(see §5.3). 

2. Initially, all the high energy ro-vibrational levels 
are considered to be in LTE. 

3. We vary T rot for all the high- J levels in each vibra- 
tional state to fit the observed line profiles. T rot is 
considered to depend on J as a 2 nd degree polyno- 
mial. 

4. In the next iteration, step 3 is repeated for the next 
ro-vibrational transition of the band under consid- 
eration. 

At each step, all line intensities are recalculated because 
a change in the rotational temperature of one of them 
implies a variation in line opacities. As the transitions 
involve higher energy levels, the variation in the opacities 
of the fitted transitions becomes lower and lower. The 
process is repeated until convergence is reached for each 
band. The remaining lines to fit belong to ro-vibrational 
transitions involving vibrational levels with higher en- 
ergies. Running the code two or three times is usually 
enough to fit the line correctly with a CPU time of less 
than 10 s per line. Unfortunately, fitting lines of hot 
bands sometimes requires repeating the fitting process 
for lines involving lower energy vibrational levels since 
the upper level of an optically thick transition can coin- 
cide with the lower level of an optically thin one. 

As discussed above, T rot may be considered equal to 
the kinetic temperature in the innermost zones for low- 
and intermediate- J levels in lower energy vibrational 
states. However, for high- J levels in high energy vibra- 
tional states, we have found that T rot is below Tk, even 
near the photosphere. In addition, for some specific tran- 
sitions, it has been necessary to include ad-hoc rotational 
temperatures in order to match the data (see Table 3). 
It is common, for two or three important ad-hoc rota- 
tional temperatures to be found in a single vibrational 
level below J ~ 30 — 35. These anomalous rotational 
temperatures are unrelated to frequencies having larger 



TABLE 3 

Ad-hoc Rotational Temperatures of C2H2 and HCN 



C2H2 

G.S.(crJ) ^4(tt s ) vs{-k u ) 


HCN 
G.S.(<7+) v 2 (tt) 


35(900,50) 26(620,200) 20(900,200) 
36(900,60) 32(425,100) 26(900,150) 
34(425,100) 32(620,200) 


20(900,250) 25(250,35) 
24(640, 85) 



Note. — Ad-hoc rotational temperatures of C2H2 and HCN 
for the lowest vibrational levels. The numbers are shown ac- 
cording to the nomenclature ^(T ro t(R dl ),T ro t(R d 2)), where J 
is the rotational level requiring an ad-hoc T ro t- The temper- 
atures are expressed in K (see the text for an explanation for 
the radial dependence). The rest of the rotational tempera- 
tures have a smooth dependence on J. We have only consid- 
ered those rotational levels involved in ro-vibrational levels up 
to J = 35 — 40 (depending on the band). Maybe new ad-hoc 
T ro t should be included for higher-J levels. 

telluric interference, and there are no apparent instru- 
mental effects. Moreover, high-J lines are in an opti- 
cally thinner part of the atmosphere than low-J ones, 
and the latter seem to be under LTE with high accuracy 
(see Figure 9). Therefore, these effects could reflect real 
physical processes such as collisional rates, radiative se- 
lective pumping effects, or overlaps of these lines with 
other spectral features (see Fonfrfa Exposito et al. 2006). 




1 10 10 20 30 
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Fig. 9. — Rotational temperatures for some C2H2 and HCN 
levels depending on the radius (left) and rotational temperatures 
for some C2H2 and HCN vibrational levels at r =R d2 depending 
on J (right). The rotational temperature is less than or equal to 
the kinetic temperature throughout the CSE (left). Fitting the 
ro-vibrational transitions sometimes needs the introduction of ad- 
hoc rotational temperatures whose dependence usually decreases 
smoothly with J (right). 

We have found that T rot for levels below J ~ 20 — 30 
follows Tk (depending on the vibrational level). How- 
ever, for high-J lines, T rot is systematically lower than 
Tk- This fact seems to indicate that the physical condi- 
tions over the inner CSE are not appropriate to maintain 
high-J levels under LTE, as we assumed in §3.2. 

5.5. Vibrational Temperatures 

The vibrational temperatures were determined from a 
fitting of line intensities, rather than by non-LTE radia- 
tive transfer calculations, but we can explain the derived 
populations with radiative pumping. 

Some excited levels of C2H2 are connected to the 
ground state only through collisions and radiative cas- 
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cades (1/4, Raman active), while others are connected 
both by collisions and direct absorption of photons (^5, 
IR active). Near the star, collisions with H 2 and He could 
play a role in the pumping of the low energy vibrational 
levels. Moreover, selection rules for C2H2 only allow ra- 
diative transitions from gerade to ungerade vibrational 
states and vice- versa (Figure 10). 



V>3 + V 5 {TT ) 



l/ 3 + l/ 4( 7T u) 



VzWf, ) 




2i/ 4 (a g + ) 



Fig. 10. — Energy of the vibrational levels of C2H2. The allowed 
radiative transitions are shown (only the perpendicular transitions 
with AU4 = ±1 and Avs = ±1, and the most intense parallel 
ones). This scheme also represents the 13 C2H2 vibrational level 
structure. The vibrational levels for H 13 CCH are similar but they 
do not express g or u. Some infrared inactive modes of C2H2 are 
weakly infrared active for H 13 CCH (i.e., the 1/4 mode described 
by Di Lonardo et al. 2002). The thin lines are allowed transitions 
while the thick ones are observed transitions. 

Consequently, the pumping of low energy vibrational 
levels could occur through absorption of near-IR photons 
emitted by both the star and the dust. These photons 
will pump C2H2 from the ground state to the strongest 
stretching mode, ^3(0^"), and its combination bands, ^3 + 
Vi{n u ), v-i + v 5 (ir g ), and v 2 + ^4 + ^5(0-+) (the notation 
is discussed in Appendix A) followed by radiative decay 
to the bending modes (Cernicharo et al. 1999). These 
infrared transitions are strong (^3(0"+)), medium (1/3 + 
^4 (?!"„), v z + v 5 (ir g )), and very strong (v 2 + va + "5 
as indicated by Mandin et al. (2005) and Herzberg II, 
p. 290. 

Radiative pumping through 1^3 (ct+) and its combina- 
tion bands, is only effective in the innermost region. Far- 
ther out, absorption of mid-IR photons coming from the 
star and the dust pump C2H2 from the ground state to 



id lo 

3 o 




II III 


11 i| 


1: 


\\ h-HH 

\ 


+ttr| \ 


1 1 1 Trnj" 
1" 



10' 
r (R„) 



10 
r (R.) 



10 



Fig. 11. — The population, P, of the vibrational levels of 
C 2 H 2 . They are labeled as follows: P = P[G.S.(cr+)], Pi = 

P[M*g)h p 2 = P[M*«)h p s = P[2M<r$)h P ± = P[ 2 MSg)}, 
P B = P[vi + u 5 {a+)}, Pa = P^4 + 1/5 (<r« )], Pt = P[^4 + u 5 (S u )}, 
P 8 = P[2v 5 (a+)], and Pg = P[2v 5 (5 g )]. The n and <5 levels are 
split into two sublevels with opposite parities, e and / (see Ap- 
pendix A). In the Figure they have been merged into a single level 
with degeneracy 2. Only vibrational levels up to 21^5 (5 9 ) have been 
plotted although higher energy levels have been considered in the 
calculations. The main box contains the vibrational levels, while 
the upper inset focuses on the high energy ones. The right boxes 
contain C2H2 density (lower) and Tpj (upper). 



z/5, and from ^4 and V5 to their combination band, V4+V5, 
and other vibrational levels, as overtones. For r > 20 R*, 
the available near-infrared photons are strongly reduced 
due to dust absorption. However, pumping through 
^5 ( i"ti) is efficient over the whole CSE. Hence, we could 
expect vibrational temperatures, T v ib, to be out of LTE 
as the distance from the star increases. 

To estimate T v ;b, we have selected some low- J lines 
from each observed vibrational band. As the vibrational 
excitation temperature, T cv , depends on T v ;b and con- 
trols the population of the upper vibrational level with 
respect to the lower one, the vibrational temperatures are 
coupled to the abundance of the molecular species under 
consideration. Hence, an iterative procedure involving 
Tk > Tev (involving T v ;b), and the abundance are required 
to derive (also fitting the observed ro-vibrational transi- 
tions) the populations of all the low- and intermediate- J 
rotational levels within vibrational levels with low and 
moderate energy, which are under LTE. So, as we will 
shows in §5.6, the first step in this procedure is to esti- 
mate the abundances over the whole envelope. The re- 
sulting T v jb for C2H2 are shown in Table 4, having been 
defined T v ib with reference to the lower vibrational state 
(e.g., T vib (G.S.) = T vib M<y - G.S.(<r+)] for C 2 H 2 ). 
Figures 11 and 12 show T ov and the derived population, 
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TABLE 4 
Vibrational temperatures of C2H2 



Vib. 
Level 


Pi 


T1 (K) 




To (TO 




Ti (K) 


G?3 


G.S. 


Po 


2330 


0.58 


900 


0.19 


685 


1.00 




Pi 


2330 


1.66 


150 


0.43 


82 


1.00 


^5(tu) 


P-2 


2330 


0.58 


900 


0.26 


625 


1.00 


2u A {a+) 


P 3 


2330 


0.58 


900 


0.58 


400 


1.00 




Pi 


2330 


0.58 


900 


0.58 


400 


1.00 


1/4 + (cri 


) Pr, 


2330 


0.58 


900 


0.58 


400 


1.00 


1/4 + u 5 (a~ 


) Pe 


2330 


0.58 


900 


0.58 


400 


1.00 


U4 + v 5 (8 u ) P7 


500 


0.98 


100 


0.16 


80 


1.00 


2^5 (ct+) 


P S 


2330 


0.58 


900 


1.40 


125 


1.00 






2330 


0.58 


900 


0.58 


400 


1.00 



Note. — Vibrational temperatures of C2H2. T\ is the vi- 
brational temperature very near the stellar surface, T2 is at the 
inner dust formation shell (r = 5.2 R»), and T3 is at the outer 
dust formation shell (r = 21.2 R*). a\, 02 , and a$ are the 
exponents of the temperature power law (r~ Q ) in Regions I, II, 
and III, respectively. Note that the vibrational temperature of 
level V4 + z/5 (S u ) is very low even close to the stellar surface. Tjj 
has been included in the last row of the Table to allow quick 
comparisons. The vibrational temperatures refer to the lower 
vibrational level. 

Pi, of the observed vibrational levels, respectively (see 
the caption of Figure 11 for the definitions of Pi). 
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Fig. 12. — Kinetic and vibrational excitation temperatures for 
the observed transitions of C2H2. T cv are merged into two differ- 
ent groups. In spite of the similarity of the dipole moments of the 
transitions, T ev [2i> 5 (a~£~) <-> v 5 (ir u )] and T ev [2v 5 (<5 9 ) i>r {-n u )] arc 
quite different from T ev [^4 + V5 (o"„ ) <-> 1/4(71-5)], T cv [u4 + ur,(8 u ) <-> 
V4,(ir g )], and T ev [v5(-ir u ) <-> G.S.(cr^)]. All the considered vibra- 
tional levels are out of LTE over the whole CSE. Note that 2^5 
and 2^5 (<5 9 ) are out of LTE even near the star. 

The analysis of the derived vibrational temperatures 
and populations (compared with those under LTE) indi- 
cates that: 

• In the innermost envelope (Region I), the popula- 
tion of the ground state is smaller than for levels 
and because they are doubly degenerate (e and 
/ parities) while the ground state is an e vibra- 
tional level. 2k and the strong radiation field near 
the star are able to pump many molecules to high 
energy levels (see Figure 11; 2i/ 5 (S g ) for example is 
at- 2100 K). 

• Viiftg) is a metastable vibrational level (see Fig- 



ure 10). It is not radiatively connected to the 
ground state through electric dipole transitions. 
Collisions and/or radiative cascades are its main 
pumping mechanisms. 

• Almost all vibrational levels, except v 5 (t: u ) and 
2v 5 (ct+), can be considered to be in LTE in Re- 
gion I. v 5 (t: u ) is out of LTE because of the strength 
of the radiative connection to the ground state 
(the dipole moment of the transition is fj, ~ 
0.31 D, Jacquemart et al. 2001). The population 
of 2f 5 (a+) departs from LTE over the whole CSE, 
even near the stellar surface (Ps/Pt- < -Ps/P-Ilte)' 
because it is connected radiatively with V5(ir u ) (the 
dipole moment of the transition 2^5 (ct+) <-> v^i^u) 
is similar to that of v^'Ku) <-> G.S.(<7+)). 

• The radial dependence of the population of all the 
vibrational levels changes in Region II with respect 
to Region I. The slope of Po and Pi falls while it in- 
creases for the other levels. The absorption of 3 jim 
radiation, which pumps C2H2 from the ground to 
high energy levels, changes the slope of the popula- 
tion of the ground state with respect to that under 
LTE throughout Regions I and II. This radiation 
also affects the v^{-K g ) level which starts to be ef- 
ficiently pumped from the radiative cascade from 
the stretching modes. The ratio P1/P0 is larger in 
Regions II and III than under LTE. However, as 
the transition V5{tt u ) <-> G.S.(<7+) is strong, the ra- 
tio P2/P1 is lower than expected under LTE. All 
the vibrational levels above ^5 are underpopulated 
in Region II. Collisions are less efficient as the vol- 
ume density in this region is low and all vibrational 
pumping is radiatively dominated. 

• 21/4(17+) and 2v i (5 g ) are metastable states like 
v/^Kg). On the other hand, v^iit^) is strongly ra- 
diatively connected to the ground state. Conse- 
quently, P3/P2 is larger than it would be under 
LTE for r >R d i. 

• Pq/P^ is larger in Region II than expected un- 
der LTE. The populations of the vibrational lev- 
els Vi + ^5(0"+), va + ^5(c~), and v 4 + v 5 (5 u ) are 
similarly increased due to radiative cascades from 
higher energy levels. However, the strong radia- 
tive transition z/4 + v^(a^) G.S.(cr+) (Hcrzberg 
II, p. 290) increases the ratio Pq/P^, due to stimu- 
lated de-excitation of molecules from i/ 4 + ^5(0-+) 
by 7.5 (im radiation. This mechanism is fur- 
ther strengthened due to the forbidden transition 
i/ 4 + i*(o--)«-G.S.(a+). 

• P4/P3, Phi Pa and P7/P6 are consistent with LTE 
as these levels are separated by an energy equiv- 
alent to 4.5, 136.1 and 10.9 K, respectively. As 
levels + ^5(0^) and i/ 4 + v^,{5 u ) are not radia- 
tively connected to the ground state, the collisions 
keep P7 ~ 2P 6 (see Figure 11). There may be a 
slight departure from LTE for the ratio P5/P4 in 
Region III (T K (100 R«) ~ 85 K), though it cannot 
be derived from the observational data. 
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• For distances larger than 20 R* (Region III), only 
the lower energy levels are significantly populated. 
The most intense transition in the outer envelope 
is V5(tt u ) <-> G.S.(<7+). Further away from the star, 
it is not possible to find any emission from C2H2. 

• As for Region II, Pi/Pq in Region III is larger than 
it would be under LTE, but it decreases rapidly, as 
in the LTE case. This comes about from absorp- 
tion of mid-infrared radiation in Region III. This 
radiation excites molecules from the state va to- 
wards va + z/5 and, then, the molecule de-excites to 
the ground state. The 3 /im pumping mechanism, 
important in Region II, loses its effectiveness at 
distances larger than Rd2 due to the decreasing in- 
tensity of the near-IR radiation. In addition, since 
the gas density and Tj< are low, collisions are inef- 
fective in de-exciting molecules from level v>A{ir g ). 
Furthermore, no molecules are in that state in the 
outer CSE. 

• Pg/Ps is consistent with LTE in Region I but 
shows a departure from LTE in Regions II and 
III. 2^5(17+) is radiatively connected to the ground 
state through 2^5(0-+) V2 + va + ^(f^") 
G.S.(<7+) and 2v 5 (a+) <- is 5 (ir u ) <- G.S.(a+). On 
the contrary, 2iy^(S g ) is connected to the ground 
state via the transition 2u 5 (5 g ) <-> V5(n u ) <-> 
G.S.(<7+) but not through 2v 5 (S g ) <-> i>2 + ^4 + 
V5(5 U ) G.S.(cr+) because V2 + va + ^(^u) 
G.S.(cr+) is forbidden. Therefore, when collisional 
pumping loses its effectiveness (r > Rdi), Pg/Ps 
becomes lower than under LTE. 

The estimation of the populations in the P3 and Pa 
levels are subject to large uncertainties because the ob- 
served lines arising from these levels are weak. We be- 
lieve that even for P3 and Pa the obtained vibrational 
temperatures are good estimators (see Figure 11). The 
rest of the vibrational populations are well determined 
based on the observational data. 

5.6. Abundances 

The condensation of refractory molecules on dust 
grains over the two acceleration zones modifies the abun- 
dances of several molecular species. It is unclear to what 
degree C2H2 (and isotopologues) condenses onto dust 
grains but our data provide a tool to check whether or 
not significant changes occur from one region to another. 

As we will show (see §8), the abundances in Region 
I has large uncertainties due to their small influence on 
most observed lines compared with the rest of the CSE. 
The angular size of the innermost region is small (~ 0'.'2) 
and, despite the high values of density and Ik, the ob- 
served flux is small with respect to that coming from 
outer shells of the envelope. On the other hand, the 
efficiency of a PSF ~ 1" is significant for the angular 
size of Region II (~ 0'.'8). Loss of light at the spectro- 
graph entrance slit, which is ~ 1'.'6 wide, is still relatively 
unimportant for Region II but is a large fraction of the 
radiation emitted in Region III. In addition, 7k and the 
density are still high in Region II (the mean Ik and n for 
C 2 H 2 are ~ 560 K and 220 cm -3 , respectively). There- 



fore, most of the emission for optically thick lines (ex- 
tending farther out than Rdi) and the bulk absorption 
for thin ones (produced in Region I) are due to Region 
II. In Region III, Ik and the density are low so emission 
from this region is anticipated to be small. Furthermore, 
Region III suffers the most from light loss at the slit, 
which pertains only to the emission component of the 
thick line profiles. Consequently Region III need only 
be considered for the very thick lines (e.g. ro- vibrational 
lines of the vibrational transition v^tTu) <-> G.S. 
for which the bulk absorption is produced by this region. 
This implies that the maximum information on the C2H2 
abundance in Region III can be obtained from absorp- 
tion in the P-Cygni profiles of these thick lines. Unfortu- 
nately, their thickness precludes the possibility of getting 
accurate information on the abundances of C2H2 at an- 
gles larger than 1" (r > 50 R»). At larger radii, the 
absorption component of the lines is almost completely 
saturated. At distances larger than 300 R», there is no 
emission from C2H2 because all the molecules remain 
in the ground state (see Figure 11) and pure rotational 
transitions are forbidden due to the lack of electric dipolc 
moment. 

The C2H2 abundances derived from the best fits for 
each region are: 

( 7.5 x 10~ 6 Region I 
^c 2 H 2 (r) = < 8.0 x 10~ 5 Region II (5) 
[ 8.0 x 10~ 5 Region III 

with an assumed initial value of 8 x 10~ 5 for the CSE. 
The derived column density is ~ 1.6 x 10 19 cm~ 2 . 

Note that the abundance in Region I is an order of 
magnitude lower than in the rest of the CSE. A higher 
abundance in that region implies diminishing its size, in- 
compatible with previous results (see §4). We will discuss 
this topic further in §8. The abundances of C2H2 in Re- 
gions II and III are similar to those obtained in previous 
studies 8 x 10 ~ 5 and 5 x 10 ~ 5 from Keady & Ridgway 
(1993) and Cernicharo et al. (1999) respectively. 

In addition, the increment in the C2H2 abundance be- 
tween Regions I and II is compatible with that yielded 
by LTE calculations (M. Agundez, private communica- 
tion) and with the results of the chemistry in a pulsating 
star published by Cherchncff (2006). The variation in 
the LTE abundance profile between both Regions ranges 
from 15% to 99% (with respect to the abundance in Re- 
gion II) , depending on the physical conditions of the very 
inner envelope (Agundez & Cernicharo 2006). The ac- 
tion of shocks on chemical abundances at radii smaller 
than 5 R* produces abundance profiles which increase 
in 99% the abundance in Region II (Cherchneff 2006). 
Both models could explain the increment of 90% derived 
in this work. We will also discuss this topic in §8. 

The uncertainties on the abundances (see Table 7) are 
compatible with a small condensation of C2H2 (and iso- 
topologues) on dust grains and the consequent existence 
of an acceleration regime between Regions II and III. 
Unfortunately, we cannot assure the existence of these 
processes with the data derived. 

6. H 13 CCH & la C 2 H 2 MODELING 

As in the case of C2H2, we can derive the abundance, 
vibrational temperatures and constrain the parameters 
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related to the physical conditions of the CSE, through 
fitting as many H 13 CCH lines as possible. The fitting 
procedure is the same as for C 2 H 2 (see §5). We have 
used the physical conditions determined previously with 
C2H2 as initial inputs for the fits of the H 13 CCH lines. 

6.1. Abundances 

The [C 2 H 2 ]/[H 13 CCH] ratio is easily determined taking 
into account that the abundance of H 13 CCH must be 
proportional to that of C2H2: 

[° 2H2 ] I 20 5 (6) 

[H 13 CCH] " 0.049 " [ ' 

where the fitted lines are shown in Table 2. This value 
is very similar to [C 2 H 2 ]/[H 13 CCH]~ 22, proposed by 
Cernicharo et al. (1996b, 1999, 2000). Hence, the abun- 
dance ratio [ 12 C]/[ 13 C] derived from our data is ~ 41 in 
IRC+10216. 

6.2. Vibrational Temperatures 

As shown in Appendix A and B, the spectrum of 
H 13 CCH presents some differences with respect to C 2 H 2 
(e.g., j/ 4 (7r)-G.S.(a+) is IR active for H 13 CCH while it 
remains forbidden for C 2 H 2 ). Consequently, although 
H 13 CCH is very similar to C 2 H 2 , its vibrational temper- 
atures are different (see Table 5 and Figure 13). The 
observed H 13 CCH transitions also support the near-IR 
pumping mechanism proposed in §5.5. For H 13 CCH, 
more transitions are allowed than for C 2 H 2 due to 
combined levels involving both stretching and bending 
modes. These differences can be understood in terms of 
the radiative pumping paths for both isotopologues (see 
below). 

• In Region I, the population of 2f 5 (cr + ) is lower than 
in LTE compared with that of ^4 + ^5(5). Radiative 
cascades from higher energy levels populate these 
states in a different manner than in the C 2 H 2 case 
(e.g., V3 + V5 <-> 2^5 is allowed for H 13 CCH), keep- 
ing level 2^5 (cr+) almost in LTE next to the pho- 
tosphere (see §5.5). Therefore, all the vibrational 
levels can be considered in LTE close to the star 
(see Figure 13). 

• The allowed transitions for H 13 CCH connecting 2^5 
to lower energy levels and the increase of the pop- 
ulation of v§ through ^4 + ^5 <-> 1/5, result in a 
decrease of T cv {2v 5 <-> i/ 5 ] in Regions II and III, 
with respect to the C 2 H 2 case. 

• In Regions II and III, P 2 /Pi, P3/P2, and P 8 /P 7 are 
lower than in LTE, as for the C 2 H 2 case. However, 
several radiative transitions in H 13 CCH, which are 
forbidden for C 2 H 2 , keep these ratios closer to LTE. 

• The v^Tig) vibrational level in C 2 H 2 is only col- 
lisionally connected to the ground state, while it 
is connected both collisionally and radiatively in 
H 13 CCH. This fact reduces the ratio P 1 /P and 
thcrmalizes the v 4 state over almost the whole CSE. 
In addition, P2/P1 is in LTE in Region I, where 
collisions can thermalize the state, but it is out 
of LTE in Regions II and III once 7k and the den- 
sity have decreased and radiative excitation has be- 
come the main pumping mechanism. As the tran- 
sition ^4 <-> G.S. is weaker than 1/5 ^> G.S. (Di 



TABLE 5 

Vibrational temperatures of H 13 CCH 



Vib. 
Level 


Pi Ti (K) 


ai 


T 2 (K) 


Ct2 j 


r 3 (K) 


a 3 


G.S. 


Po 


2330 


0.58 


900 


0.56 


410 


1.00 


1/4 (71- ) 


Pi 


2330 


0.58 


900 


1.92 


60 


1.00 




P 2 


2330 


0.58 


900 


1.43 


120 


1.00 


2v 4 (cr + ) 


Pi 


2330 


0.58 


900 


0.58 


400 


1.00 


2u 4 (S) 


Pi 


2330 


0.58 


900 


0.58 


400 


1.00 


1/4 + 1/5 (cr 4 


■)fls 


2330 


0.58 


900 


0.58 


400 


1.00 


U4 + 1/5 (a~ 


)Pe 


2330 


0.58 


900 


0.58 


400 


1.00 


v 4 + u 5 (8) P 7 


2330 


0.93 


500 


2.78 


10 


1.00 


2v 5 (<r+) 


Ps 


2330 


0.58 


900 


0.58 


400 


1.00 


T K 




2330 


0.58 


900 


0.58 


400 


1.00 



Note. — Vibrational temperatures of H 13 CCH. See Ta- 
ble 4 for details about the meaning of the constants. 

Lonardo et al. 2002), we can conclude that pump- 
ing of molecules from the ground state to v 4 and 
1^5 through radiative cascades from high-energy vi- 
brational levels seems to be almost the same (see 
Figure 10). 

• P 3 j Pi is lower than in LTE in Regions II and III for 

several reasons: 2^4 (er + ) is radiatively connected 
to ^4(71") and G.S. (cr + ); molecules are pumped 
from the ground state to V5(ir) through radiative 
cascades, and collisional pumping is not effective 
enough to keep the vibrational populations under 
LTE in these regions. 

• As for C 2 H 2 , the model points out that P4/P3, 
P5/P4, and P7/P6 can be considered to be in LTE. 
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Fig. 13. — Kinetic and vibrational excitation temperatures for 
the H 13 CCH observed transitions. As for C2H2, T ov are merged 
into two different groups, although in this case, the transition 
^* G.S. (cr+) stands alone with a T ev not much lower than 
Tjj. All the selected transitions are under LTE in Region I. 



T ev [^ 5 (7r) <-> G.S.(cr + )] < Tk because the transi- 
tion is strong. On the other hand, the transition 
^4+^5 ^5 could depopulate which is com- 

patible with Tcv[^4 + ^5 ^4] < T cv [iy 5 ^ G.S.], in 
contrast to the behavior of these temperatures for 
C 2 H 2 where they are equal. 
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6.3. Rotational Temperatures 

Unlike C2H2, all H 13 CCH chosen lines can be fitted 
quite well with the rotational levels in LTE. No ad-hoc 
rotational temperatures are needed. 

6.4. Searching for 13 C2H2 

Given the large abundance of C2H2, detection of 
13 C2H2 ro-vibrational lines may be possible. The dou- 
ble isotopic substitution of the C-atoms greatly reduces 
the abundance of this species compared with C2H2. As 
we have derived in §6, [ 12 C]/[ 13 C] ~ 41. Consequently, 
[ 13 C 2 H 2 ]/[C 2 H2] ~ (1/41) 2 - 6 x 10- 4 . Unfortunately, 
the strongest bands of 13 C2H2 are in a range of the spec- 
trum crowded with strong C 2 H 2 , H 13 CCH, and HCN 
transitions (see Figure 1). We have compared the optical 
depth of one of the strongest 13 C 2 H 2 lines, v 5 (ir u )R e (3) , 
with the same line for C2H2. The dipole moment em- 
ployed in the calculations for 13 C2H2 has been set equal 
to that of C2H2 due to the lack of data on intensities 
quoted in the literature. The optical depth of this line 
has been computed by running the code with the param- 
eters derived from the fits to the observed C2H2 lines and 
using the double-substitution isotopic ratio given above. 
The derived optical depths along the line of sight, with- 
out considering the Doppler effect, are 630 for C2H2 and 
0.033 for 13 C 2 H 2 . The low optical depth for the 13 C 2 H 2 
line results in a maximum absorption < 1% of the contin- 
uum flux and cannot be observed in our data. An equal 
or smaller absorption is expected for other lines of the 
same band. 

7. HCN AND H 13 CN, MODELING 

Modeling HCN (and H 13 CN) lines adds additional con- 
straints to the physical conditions over the CSE for 
r < 100 R*. The spectra shown in Figures 1—6 indi- 
cate the presence of many ro-vibrational lines from the 
v 2 bending mode and several associated overtones. As for 
C2H2 and H 13 CCH, we have selected a sample of HCN 
(and H 13 CN) lines to fit to the data and to derive phys- 
ical and chemical conditions through the CSE (see some 
fits in Figure 14). The observation of transitions arising 
from overtones and H 13 CN allows us to study optically 
thinner lines compared to those of the fundamental bend- 
ing mode. 

The differences between H 13 CN and the main isotopo- 
logue, HCN, lead to a significant increment in the molec- 
ular mass and a slight variation in the interatomic dis- 
tances. Consequently, the wavefunctions of each molec- 
ular state and their related energies are not equal for 
HCN and H 13 CN. Although some spectroscopic work on 
H 13 CN can be found in the literature, there is no data 
available about the intensities of its transitions, while 
HCN has been studied more extensively. Therefore, we 
have used the dipole moments of HCN for H 13 CN as a 
reasonable approximation for the transition strengths. 

The number of lines used for the fits is shown in Ta- 
ble 2. Information about frequencies, intensities, and 
nomenclature for HCN and its isotopologucs can be 
found in Appendixes A and B, and references therein. 

7.1. Vibrational Temperatures 

As seen in previous sections, C2H2 does not present any 
permanent dipole moment (it is very small for H 13 CCH 



TABLE 6 

Vibrational temperatures of HCN and H 13 CN 



Vib. 
Level 


Pi 


Ti (K) 




T2 (K) 


a 2 '. 


f3 (K) 


C*3 


G.S. 


Po 


2330 


0.58 


900 


1.56 


130 


1.00 


V2 (-71- ) 


Pi 


1250 


0.82 


600 


0.62 


250 


1.00 


2v 2 (o+ 


)P2 


800 


1.26 


100 


1.24 


17.5 


1.00 


T K 




2330 


0.58 


900 


0.58 


400 


1.00 



Note. — Vibrational temperatures of HCN and 
H 13 CN. See Table 4 for details about the meaning of 
the constants. The parameters of the vibrational tem- 
peratures of H 13 CN are only those related to the ground 
state. 



and can be neglected to a first approximation). Consid- 
ering rotational LTE as suitable initial input, there will 
not be any pure rotational emission. Fortunately, in spite 
of the permanent dipole moment, HCN (and H 13 CN) 
seems to be under LTE for low- J rotational levels, mak- 
ing easier the determination of the abundance and T v ib 
(see §7.2). 

Near R d i, T d ~ 800 - 850 K and the dust black-body 
emission peaks at ~ 3.5 ^m (~ 2850 cm -1 ). This radia- 
tion pumps HCN from the ground state to the stretching 
and combination levels v^{a + ), v 2 + ^3(71"), and ^i(cr + ), 
with energies 2096.85, 2807.05, and 3311.48 cm" 1 respec- 
tively (see Figure 15). In addition, several rotational lev- 
els of 4^2 (ct + ) (with energy 2802.96 cm -1 ) are strongly 
connected to the ground state through a resonance with 
the corresponding rotational levels of v 2 + 1/3(71"). This 
resonance could explain the HCN masers observed by 
Schilke & Menten (2003) in the Av 2 ((j + ) state. Radiative 
deexcitation changes the population of mode v 2 (n) and 
its overtones with respect to the LTE case. 

The vibrational temperatures derived from the fits are 
presented in Table 6 and plotted in Figure 16. As in 
the case of C2H 2 , we have labeled the populations of 
the levels as follows: P = P[G.S.(ct+)], Pi = P^tt)], 
P 2 = P[2z/ 2 (a+)], and P 3 = P[2v 2 {5)\. 

• T v - 1 \ } [v 2 {'k) <-> G.S. (cr + )] ~ Tk in Region I, depart- 
ing from Tk in Regions II and III. In Region I, col- 
lisions thermalize ^2(71"), while at larger distances 
only infrared photons play a role in the pump- 
ing. Vibrational levels tr(o" + ) and 1/3 (cr + ) are ra- 
diatively connected to the bending level v 2 (n). In 
addition, v 2 {ii) can also be populated through the 
loop G.S. (er + ) — > v 2 + ^3(71") — > v?,{o- + ) — > v 2 {tx). 
The band v 2 (k)~ G.S. (cr + ) is quite strong and the 
deexcitation rate is high. 

• 2v 2 (o- + ) and 2u 2 {8) are out of LTE over the whole 
CSE, as well as close to the stellar photosphere. 
The dipole moment of the transitions involving 
these levels are so high that they are not ther- 
malized even in the dense warm innermost region. 
Changing T v ;b of these levels at the stellar sur- 
face from 1250 K to 2330 K produces a change of 
20 — 30% in the emission of 2v 2 [a Jr ) <-> v 2 (-n:) ro- 
vibrational transitions. Hence, our data is very sen- 
sitive to the value of T v ib. 2^ 2 (cr + ) can be pumped 
from the ground state and from v 2 + v 3 {a + ). How- 
ever, the 2v 2 (5) level is not connected to the ground 
state except through the radiative cascades from 
higher energy levels (see Cernicharo et al. 1999). 



The Dust Formation Zone of IRC+10216 



21 



dv? —v? R„ 



T 3 l-O 




v(km s 1 ) 

Fig. 14. — Fits of some lines of HCN. The black lines are the observed spectra while the red lines are the fits. The data have been 
corrected from observational frequency deviations adding/suppresing the mean difference between the observational and the laboratory 
frequency of many lines (see §2). The existing discrepancies are small and can be assumed as blends with other lines and/or observational 
errors (the spectroscopic notation is described in Appendix A). 



Therefore, P3/2P2 is always lower than 2P2/P1 (we 
add the factor 2 to account for the difference in the 
degeneracy between levels 2z/ 2 (cr + ) and 2u 2 {8), and 
2u 2 (a + ) and 1*2(71-)). 

T v ih[v2— G.S.] and Tvit,[2i/ 2 — V2] are equal to each 
other at ~ 10 R* (see Figure 16). Outside of the 
region where the NIR radiation emitted by the dust 
arises, the pumping via higher energy levels from 
the ground state (through 7 /jm radiation, which 
is the wavelength corresponding to the maximum 
emission of a black-body at Tbb — 400 K) be- 
comes less important, meaning that the population 
of level 2v 2 {cr + ) increases with respect to 1*2 (7r). 



7.2. Rotational Temperatures 

The rotational structure of the ground state seems to 
be in LTE for rotational levels up to J = 20. It is nec- 
essary to adopt ad-hoc rotational temperatures for rota- 
tional levels J = 20 and 24, when fitting the ^ 2 R e lines 
(see Table 3). The line z/ 2 (7r)-G.S.(o-+)R e (21) can be 
fitted more accurately with an ad-hoc T rot at r = Rd2 



for J = 20 (although the fit is relatively good even un- 
der LTE). However, this ad-hoc temperature is neces- 
sary to fit higher- J ro-vibrational lines. J = 24 must be 
out of LTE. Under LTE, the absorption and the emis- 
sion of the line ^ 2 (vr)-G.S.(cr+)R e (25) is 15% and 70% 
larger respectively compared with the observed feature. 
In addition, considering that this level departs from LTE 
with the same smooth dependence as adjacent levels, 
does not give a good fit as synthetic profile presents 77% 
more emission and 6% more absorption than the best 
fit. Using an ad-hoc rotational temperature improves 
the fits to higher- J ro-vibrational lines. The fits suggest 
a quadratic deviation from LTE for the rotational pattern 
of the ground state beginning at J = 20 with a decrease 
per rotational level of 96 — 7.8( J — 20) K at r = R c n 
and of 46 - 3.2( J - 20) K at r = R d2 . From the fits of 
bands 1*2(71") — G.S.(ct + ) and 2z/ 2 (c T+ ) — ^(tt), we can as- 
sume that 1/2 (tt) is in LTE for the rotational levels up to 
J = 22 departing smoothly from LTE for higher levels. 
The most important ad-hoc temperature corresponds to 
J = 25, being necessary to fit ro-vibrational lines involv- 
ing higher- J levels. We can assume that 2i/ 2 (a + ) rota- 
tional levels are in LTE for J < 26, with it necessary to 
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Fig. 15. — Energy of the vibrational levels of HCN. This pattern 
is the same for the H 13 CN isotopologue. The transitions given 
are the perpendicular ones with Ai>2 = ±1 and the most intense 
parallel ones. The thin lines are allowed transitions. The thick 
lines are observed vibrational transitions. The red dashed double- 
arrow indicates the resonance between several rotational levels of 
the vibrational states + v%(t{) and Au2(u + ). 
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Fig. 16. — Kinetic and vibrational excitation temperatures for 
HCN observed transitions. See text for details. 

add an ad-hoc rotational temperature for J = 27. Nev- 
ertheless, fitting the bands Iv^iS) — v?,(tk) also requires 
anomalous non-LTE rotational temperatures at the pho- 
tosphere Rdi, and Rd2 for lines J = 12 and 26. The 
best fits to all the lines involving anomalous rotational 
temperatures require T rot (Rd2) < 50 K, although tem- 
peratures below 100 K are acceptable despite poorer fits. 



7.3. Abundances 

Low-J HCN lines of band ^(t)— G.S((T + ) are opti- 
cally thick due to the large dipole moment of the vi- 
brational transition. Many lines are saturated and little 
information can be obtained from them. In particular, 
the abundance is a parameter which can not be accu- 
rately determined through fits of low energy HCN lines. 
The intensity of the G.S.(cr + ) pure rotational lines is di- 
rectly related to the abundance of HCN while the inten- 
sity of ro- vibrational lines also depend on T V ib, which is 
initially unknown. The lack of radio observations of these 
Regions hinders an accurate determination of the HCN 
abundance. Hence it is better to obtain the H 13 CN abun- 
dance by fitting the optically thin ^(tt) — G.S.(cr + ) lines, 
and using the derived ratio [ 12 C]/[ 13 C]~ 41 to determine 
the abundance of HCN. We found that the abundances 
which produces the best fits are: 

( 1.23 x 10~ 5 Region I 
ihcn(0 = < 4.5 x 10~ 5 Region II (7) 
[ 2.0 x 10~ 5 Region III 

which means that the column density is 1.6 x 10 19 cm~ 2 , 
i.e., similar to that of C2H2. In fact, HCN seems to be 
more abundant than C2H2 in Region I, where the gas 
density is rather high. However, as in the C2H2 case, the 
inner abundance is poorly determined (see §8) and the 
middle and outer ones are compatible with those found 
in previous works (Keady & Ridgway 1993; Cernicharo 
et al. 1999, 1996b; Wiedemann et al. 1991; Dayal & Bieg- 
ing 1995; Lindqvist et al. 2000). The abundance of HCN 
between Regions I and II increases by a factor 3.7 while 
for C2H2 the increase is by a factor of 10, incompatible 
in this case with the predictions of chemical LTE mod- 
els (M. Agundez, private communication) and non-LTE 
stellar pulsating ones (Cherchneff 2006), where a dimin- 
ishing in the HCN abundance with growing radii is found 
(see §8 for a discussion on this topic). The decrease by a 
factor of ~ 2 in the abundance between Regions II and 
III seems to be real (see Table 7). 

8. SENSITIVITY TO DIVERSE PARAMETERS AND 
DISCUSSION 

To obtain complete information about the CSE, it is 
necessary to study the behavior of the model while vary- 
ing the parameters and their uncertainties. For this pur- 
pose, we have selected lines of all the molecules studied 
in the present work. Each line is radiatively active over 
different regions of the CSE. For example, the H 13 CCH 
line 21/5(6) <-> ^5(7r)R e (7) extends until 22 R* and has an 
optical depth of 8.50 x 10 -2 in Region I, 1.41 in Region 
II, and 9.94 x 10~ 3 in Region III. On the other hand, 
the C2H2 line 1^5(7^) <-> G.S.(tr+)R e (30) extends until 
27.5 R* and has an optical depth of 3.39 in Region I, 
39.0 in Region II, and 1.35 in Region III. The first line 
is useful to study Region II while the second is useful for 
Region III because they are optically thin. Table 7 shows 
the calculated errors and the lines used for this purpose. 
The chosen lines appear to be the best for a determina- 
tion of the errors in the physical parameters. The veloc- 
ity gradients introduce different spatial contributions to 
the line profiles. Hence the impact of each region of the 
CSE becomes measurable through the analysis of these 
lines. Moreover, a line that is globally optically thick 
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can also carry information on specific regions of the en- 
velope. However, there are some parameters such as Tdi 
or Rdi whose errors must be calculated by fitting only 
the continuum because the changes they produce on the 
lines, when varied, can be overcome by modifying other 
parameters, e.g., vibrational temperatures. 

The results obtained in this paper have been based 
on several hypotheses, most of them supported by 
observational data. However, the assumption for 
which the observations of IRC+10216 by ISO/SWS and 
IRTF/TEXES are compatible (see §4) needs discussion. 
As we pointed out in §2, the difference in the IR phase be- 
tween both observations is Aipm ~ 0.3. Hence, the phys- 
ical properties in the innermost CSE are different for ob- 
servations of gas (IRTF/TEXES) and dust (ISO/SWS). 
Each pulsation of the star is followed by an increment 
in the emitted radiation and, consequently, Rdi and Tdi 
are magnified because of the new input of optical and IR 
photons from the star and the ejected gas. Nevertheless, 
the large optical depth of dust shells near Rdi diminishes 
the number of short wavelength photons from the stel- 
lar surface as the radius increases. The rest of the high 
frequency stellar emission is largely diluted and docs not 
significantly affect the middle and outer envelope. The 
result is that only the inner shells of the dusty CSE are 
significantly heated. The dust grains near R d i reemit 
the absorbed energy at longer wavelengths increasing 
the observed NIR intensity of the continuum. Since FIR 
and radio emission from the source arise almost entirely 
from the outer dust shells (the black-body emission at 
Tbb = 150 K peaks at 19.3 ^m and the dust shell at 
that temperature at ~ 450 R*), stellar pulsation has lit- 
tle effect on the long wavelength range of the continuum. 
However, it does affect the MIR contribution (Monnier 
et al. 1998). By modifying T dl from 850 to 950 K and the 
exponent ay from 0.39 to 0.42 to model the increment 
of Td in the inner dusty CSE during a stellar pulsation, 
the continuum changes less than 15% over the frequency 
range considered in this paper. Any variations intro- 
duced in the lineshapes can be overcome by changing the 
T ov profiles. Despite these changes, the rough behavior 
of T cv remains. The abundance profile does not seem to 
be affected by these modifications of Td. The increase in 
the abundance detected between Regions I and II seems 
to be real and not an effect of an incompatibility between 
the ISO/SWS and IRTF/TEXES observations. 

The expansion velocity in Region I derived by us is 
larger than that suggested by Keady et al. (1988) and 
Keady & Ridgway (1993) by a factor of 1.75. The main 
effect produced by this difference on the P-Cygni pro- 
files is found at the red wing of the line emission, at 
positive velocities. Increasing the expansion velocity of 
the gas in Region I expands the red wing of P-Cygni pro- 
files. A reduction to 3 km s _1 does not improve the fits. 
Expansion velocities larger than 5 km s _1 seem to give 
better results, but are improbable. The synthetic profiles 
with 5 km s _1 fit the observed lines quite well, although 
a slight lack of emission at velocities near terminal in 
C 2 H 2 z/ 5 (7r u )-G.S.(a+)R e and HCN ^ 2 (7r)-G.S.(cr + )R e 
lines can be seen. 

In order to detect any effect produced by the differ- 
ence between the pulsation phase of the IRTF/TEXES 
and ISO/SWS observations, we have varied the param- 



eters Rdi and Td. They do not cause any modification 
of the width of the red wing. The fact that the width 
of the red part of the observed P-Cygni profiles is larger 
than that of the synthetic ones, considering the accepted 
microturbulence, i.e., Avi — 5 km s _1 , suggests either 
more emission from the outer CSE only at the back of the 
star or larger linewidths in the innermost envelope. The 
former explanation is in disagreement with the spheri- 
cal symmetry established by a large number of observa- 
tions. The latter seems to be more realistic. A possi- 
ble explanation could be an onion-structured Region I 
with shells expanding at positive and negative velocities 
corresponding to expansion and collapse (Bowen 1988). 
However, we could suggest another scenario explaining 
large linewidths in the innermost CSE where it would be 
a clumpy region. The clumps would be composed of hot 
gas and would move inwards and outwards at high veloci- 
ties along different radial directions. Unfortunately, little 
information on the spatial distribution of these clumps in 
Region I could be derived through the observed spectrum 
due to heavy dilution. Both scenarios have been invoked 
by Fonfrfa Exposito et al. (2006) in order to explain the 
observed SiS masers towards IRC+10216, where each 
maser feature has widths of ~ 2 — 3 km s _1 provid- 
ing valuable information about the physical conditions 
of the emitting regions, i.e., shells or clumps. Modeling 
these structures as in §5.2 with Av\ = 30 km s _1 and 
I = 1.5 R* produces synthetic profiles having an emis- 
sion wing wider than with Av\ — 5 km s -1 . Related to 
the derived column densities for C 2 H 2 and HCN, they 
are beam averaged and thus their actual values for the 
clumps in Region I could be much larger if we take into 
account dilution. More spatial resolution is needed to 
study the complex, small scale structure of the inner- 
most CSE. 

The possible condensation of C 2 H 2 (and isotopologues) 
onto the dust grains in the outer acceleration zone is 
somewhat constrained by our derived uncertainties in 
the abundances. In Regions II and III, these allow for 
20 — 30% of the C 2 H 2 to be deposited onto grains. We 
can estimate how much the dust grains should increase 
in radius from this deposition. Using the density of C 2 H 2 
and the density of dust grains at Rd 2 derived from the 
fits (~ 1200 cm~ 3 and ~ 3 x 10~ 4 cm~ 3 respectively), 
assuming that 30% of the C 2 H 2 turns into solid state 
with density ~ 1 g cm~ 3 and that the dust grains are 
spheres of radius 0.05 /xm, we calculate that 1.2 x 10 6 
molecules of C 2 H 2 condense onto a single dust grain. 
The contribution to the mass of the dust grain is 10% 
and the radius of the grain grows by 3.2%. For a larger 
initial dust grain radius, the increase is even smaller. The 
small increase in the size of the dust grains suggests that 
the condensation of C 2 H 2 onto the dust grains does not 
contribute significantly to the grain growth and so other 
carbon molecules are more likely to be responsible for 
it. However, it might be possible that C 2 H 2 reacts with 
other molecules on the grain surface giving rise to more 
complex organic molecules. 

Regarding the decrease of the HCN abundance in Rd 2 , 
it might be thought that HCN suffers a depletion onto 
dust grains. In that case, it would be reasonable to ex- 
pect some emission/absorption in the observed contin- 
uum at wavelengths between 4 and 5 (im due to the 
stretching mode of the group C=N (see, for example, 
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TABLE 7 

Errors 



Abundance 



Region 



Value 



Parameter 



Unit 



Value 



x(C 2 H 2 ) 



x(H 13 CCH) 



x(HCN) 



z(H 13 CN) 



I 

11 
III 
I 

11 
III 
I 

11 
III 
I 

II 
III 



( 




(5.0t?:°*) 



x 10~ 7 



Rdi 

R d2 

Tk(Rdi) 
Tk(R d2 ) 
"i 

T vib (i/ B R e ,R dl ) 

7 vib( I/ 5R.e,Rd2) 

r d i 

"e,I 

«e,II 

f e ,III 

a d 
r(ll /im) 



R* 
R* 
K 
K 



K 
K 
K 

km s -1 
km s -1 
km s -1 



+0.6i 
0.5 



5.2 
21 ±3 b 



900 



+300 c 
-200 

+27b 



400_ 25 
0.58 ±0.16° 
0.58^ 
175 1 



+0.05 b 
0.04 

+22 b 
20 

82 +i0b 



'-9 

850 ± 50 



5.0 



+1.6 e 
0.9 



n.otij 

14.5 ± 0.4 C 
0.39 4 



0.70 



+ 0.08 i 
0.06 
+0.13i 
0.11 



Note. — The parameters derived by the model with their errors, estimated through a sensitivity study. All the errors 
translate into an interval of the normalized flux with a maximum width of 20%, containing the best fit. The uncertainty 
of a given parameter has been calculated by comparing the synthetic spectrum with the observed one and by forcing the 
former to be within that interval, while the rest of the parameters are modified until getting the maximum/minimum value 
for the considered one (see §3.2 for more information about the method). Some of the values in the upper Table could be 
slightly different from those appearing along the text, since the latter produce the best fits and the former produce just 
acceptable fits and are complemented with the uncertainties of the parameters. The x's are the abundances in the three 
regions of the envelope. R^i an d R-d2 are the positions of the inner and outer dust formation shells, respectively. The kinetic 
temperature, Tj£, has been evaluated at R^i and R^- a l an d an are the Tj^ exponents in Regions I and II. T v ;b,(^5Re,Rdi) 
and T v ;{ 5 (^5R e ,R ( j2) a rc the vibrational temperatures of the transition 1/5 — G.S.R e evaluated at R dl an d Rd2 respectively 
(shown in the Table as an example of the T V ; D uncertainties). T^i is the dust temperature at R^i- The expansion velocity 
of the gas, v e , is shown over the three regions of the CSE. ay is the exponent of the dust temperature law. Finally, the last 
parameter is the dust optical depth, r, at 11 fim. 

A list of the lines used to calculate the error for each parameter is shown below. Each line is accompanied by the set of 
parameters P = [r x % (R*), Total > T' HI' T IIl]> where r x % is the minimal radius at which the opacity of the line is 1% of the 
maximum, Total is the total optical depth, and tj, tjj, and tjjj are the optical depths in Regions I, II, and III, respectively. 

a C 2 H 2 2u 5 (8 g ) ^ 5 (7r„)R e (30) with P = [20.1,3.57,0.549,2.83,0.00664] 
b C 2 H 2 u 5 (tt u ) <-► G.S.(cr^)R e (30) with P = [28.4, 42.0, 1.83, 37.2, 0.988] 
C C 2 H 2 u 5 (tt u ) <-> G.S.(cr+)R e (6) with P = [206, 241, 1.60, 103.0, 127.4] 
d H 13 CCH 2z/ B (5) <-> ^ 5 (7r)R e (7) with P = [22.6, 1.47, 0.0492, 1.38, 0.00660] 
e H 13 CN v 2 (-k) <-> G.S.(<7+)Re(14) with P = [46.4, 3.14, 0.132, 2.64, 0.234] 
f H 13 CCH u 5 (tt) <-> G.S.(<r+)R e (26) with P = [37.8,2.98,0.0629,2.59,0.188] 
SHCN u 2 {tt) <-> G.S.(<r+)R e (22) with P = [32.8,62.6,5.16,53.9,1.07] 

h The absolute error of the abundance in the given Region for the isotopologue considered has been calculated with the 
relative error for the other isotopologue. Taking the ratio [ 12 C]/[ 13 C] as constant throughout the CSE the relative errors of 
both isotopologues must be equal. 

'The errors have been calculated by fitting the continuum with different values of the parameter. 



Lacy et al. 1984; Pendleton et al. 1999). Nevertheless, 
no significant emission/ absorption can be seen in the 
ISO/SWS spectrum of IRC+10216 at these wavelengths 
although it is possible that some HCN molecules con- 
dense onto dust grains. However, their contribution to 
the continuum could be masked by the CO vibrational 
band at 4.67 /im. Even in that case, we think that this 
deposition process does not account for the abundance 
decay observed at Rd2, which could actually be due to 
chemical processes. The most important HCN decrease 
is produced by photodissociation in the external layers 
of the CSE (Agiindez & Cernicharo 2006, and references 
therein) . 

9. CONCLUSIONS 

IRC+10216 has been observed from 11 to 14 itm with 
the high resolution spectrograph TEXES at the 3 m 
IRTF. We complemented this data with information on 



the continuum observed with ISO/SWS. We have iden- 
tified 462 ro- vibrational lines of C2H2 (involving vibra- 
tional levels up to 3^5 at ~ 2185 cm -1 ), 95 of HCN (in- 
volving levels up to 2v 2 at ~ 1425 enr 1 ), 106 of H 13 CCH 
(involving levels up to 2v 5 at ~ 1454 cm -1 ) and 7 of 
H 13 CN (the fundamental band v 2 at ~ 707 cm -1 ). By 
means of a model of an AGB star developed by us, we 
have fitted over 300 lines. The results can be summarized 
as follows: 

1- The geometrical structure and physical properties 
of dust and gas over the whole envelope are compatible 
in many cases with those already proposed. However, 
the values of certain parameters are quite different. In 
particular, we find the position of the outer acceleration 
zone (~ 21 R» ~ 0'.'4) to be farther out than found by 
previous work. 

2.- The abundances of C2H2 and HCN in the inner- 
most CSE, Region I (8 x 10~ 6 for C 2 H 2 and 1.2 x 10~ 5 
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for HCN with [ 13 C]/[ 12 C]~ 41) are lower than those over 
the outer envelope, Regions II (8.0 x 1CP 5 for C2H2 and 
4.5 x 1CT 5 for HCN) and III (8.0 x 10~ 5 for C 2 H 2 and 
2.0 x 10~ 5 for HCN). For HCN, the ratio of the abun- 
dance in Region II to that in Region I is only a factor 
~ 4 but, in the case of C2H2, it is a factor ~ 10. The 
latter result is in accordance with the increase predicted 
by the chemical models. Contrarily, the derived abun- 
dance of HCN grows as radius increases while the chem- 
ical models suggest a decay. In addition, the determined 
abundances in Regions II and III allow the condensation 
of < 20 — 30% of molecules of C2H2 onto the dust grains. 
On the other hand, the decrease of the HCN abundance 
when reaching the outer dust formation zone is probably 
due to chemical processes. 

3. - The vibrational temperatures determined by fit- 
ting the lines suggest the existence of a complex pump- 
ing mechanism driven by near-IR radiation (3.5 — 7 /im) 
emitted by the star and the inner dusty CSE. High- 
energy vibrational levels play an important role in the 
non-LTE pumping of low- and mid-energy vibrational 
states throughout Region II and the edge of Region 
III through radiative cascades. Reliable modeling us- 
ing spectroscopic methods of any molecular species in 
the innermost CSE requires an analysis of high-energy 
vibrational levels at temperatures as high as ~ 3000 K. 

4. - Most of the rotational levels behave as expected: 
low-J rotational levels can be considered to be in LTE, 
while high- J ones are not thermalized. Interestingly, we 
have found it necessary to add ad-hoc rotational temper- 
atures to fit several ro-vibrational lines. These T rot are 
quite different from those of adjacent levels. The involved 
ro-vibrational transitions do not seem to be affected by 
instrumental or telluric effects so we assume these vari- 
ations to be real and produced by molecular processes 
such as overlaps with lines of other molecular species. 

Infrared spectroscopic data provide us with extensive 
information about physical conditions of warm sources 
with good spatial resolution and without resorting to in- 



terferometric methods. However, physical processes and 
the time-dependent chemistry in the innermost CSE re- 
main unknown. Future high angular resolution observa- 
tions with TEXES in the infrared domain will provide us 
with high quality data to study the dynamics, chemistry, 
and physical conditions in the warm asymmetric inner- 
most envelope. On the other hand, ALMA will supply 
us in several years with the interferometric observations 
needed to delve more deeply into the dynamics of the 
near environment of the central star and to further our 
understanding of dust formation and growth. 
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APPENDIX 
LINE FREQUENCIES 

The spectroscopic data for all the molecular species have been taken mainly from the HITRAN Database (Rothman 
et al. 2003). Further information came from data published by Herman et al. (2003), Kabbadj et al. (1991), and Di 
Lonardo et al. (1993, 2002) for C2H2 and its isotopologues, and from Maki et al. (1996, 2000) and Malathy Devi et al. 
(2005) for HCN and H 13 CN. The spectroscopic data relative to the modeled lines can be found in Table Al. 

The molecules studied in the present work, acetylene and hydrogen cyanide, are linear. Consequently, the notation 
we have used is the same for both species. The notation adopted to refer to the vibrational normal modes is Vih>i{£ p s ), 
where i is the mode number (i = 1,2,..., 3 A — 5 — d, where N is the number of atoms of the molecule and d is the 
number of degenerate modes), Vi is the vibrational quantum number related to the z th normal mode, £ is a Greek letter 
corresponding to the quantum number of the total vibrational angular momentum of the vibrational state (I — = a, 
£=1=tt,£ = 2 = 5, etc.), S is the symmetry of the ro-vibrational state with respect to the molecular midplane (S = g 
for gerade or even states and S = u for ungerade or odd ones), and p is the parity of the total molecular wavefunction 
(+,—). If no parity indicator appears, the label refers to the + and — levels simultaneously. If no symmetry indicator 
is present, the molecule is not symmetric with respect to the molecular midplane, as is the case of H 13 CCH and HCN. 
For the combination bands, t, S, and p refer to the mixture of all vibrational modes which participate in the combined 
state. Combination bands are sometimes split in different vibrational levels having the same quantum number £ =/= 0, 
each one with both parities + and — . In these cases, Roman numbers are added to label these states from higher to 
lower energies (see Figure 10). In addition, the interaction between the vibrational and molecular angular momenta 
splits the rotational levels of vibrational states with £ ^ into two sub-levels with opposite parities denoted e and /. 
The energy gap between those sub-levels for a given rotational state depends on J and on the vibrational state. For 
low-J levels it is between 200 and 300 times lower than the rotational constant for acetylene and hydrogen cyanide. 
The gaps are large enough to be observed in high-resolution spectra, as in our case. In the label for each transition, 
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TABLE Al 

Spectroscopic Data of the Modeled Lines 



Mol. 


Transition 


Preq. (cm x ) 


R 2 (10- 2 D 2 ) 


Ei (cm- 1 ) 




9u 


9l 


HCN 


2u 2 (<5) <-> y 2 (7r)R / (3) 


726.423637 


7.460 


729.8054 


37.500 


9 


7 


HCN 


2u 2 (a-+) <-> i/ 2 (7r)R e (8) 


726.700075 


3.740 


818.3829 


40.000 


19 


17 


H 13 CCH 


1/5 (tt) «-> G.S.(<T+)R e (l) 


732.820280 


8.744 


2.2969 


15.000 


5 


3 


H 13 CCH 




734.177200 


5.118 


671.4353 


35.000 


17 


15 


H 13 CCH 


1/5 (tt) «-> G.S.(<7+)R e (2) 


735.115630 


6.995 


6.8907 


20.000 


7 


5 


HCN 


2^(5) <-► i . 2 ( 7r )R / (6) 


735.315339 


7.480 


774.3662 


51.429 


15 


13 


C 2 H 2 


1/4 + ^5 (<4~) <-> Vi(-Kg)R e (7) 


735.543410 


5.118 


677.5067 


35.000 


51 


45 


HCN 


u 2 (tt) «-> G.S.(<r+)R e (7) 


735.611573 


3.740 


82.7713 


45.000 


17 


15 


H 13 CCH 


1/5 (tt) «-> G.S.(<7+)R e (3) 


737.410380 


6.247 


13.7813 


25.000 


9 


7 


C 2 H 2 


IA4 + ^(fj ) <-> I^4(7r S )Re(8) 


737.980170 


5.176 


696.3085 


40.000 


19 


17 



Note. — Spectroscopic data relative to the modeled lines. The transitions are labeled in Columns 1 
and 2, while in Columns 3 to 8 are shown the frequency in cm -1 , the square dipole moment in 10" 2 D 2 , 
the energy of the lower level in cm -1 and the line strength of each transition, and the degeneracy of the 
upper and lower levels respectively. See the text for details on the notation relative to the transitions. [The 
complete version of this table is in the electronic edition of the Journal. The printed edition contains only 
a sample.] 

the parity e or / of the lower ro-vibrational level is included. Deriving the parity of the upper level can be done with 
the selection rules governing transitions in linear molecules: 



Av, = ±1 
Av t = ±1,±2,±3, 
A£ = 0, ±1 

9 u 
g <-> g u u 
A J = 0,±1 

e <->■ e / <-> / 
e / 

where the transitions with Al 



i = bending modes, 

i = stretching modes, 

Ai = allowed only for £\ ov , = 0, 

forbidden, 

allowed (C 2 H 2 , 13 C 2 H 2 ), 
forbidden (C 2 H 2 , 13 C 2 H 2 ), 

A J = forbidden for J = J = 0, when Al = and 4, 

for parallel transitions; allowed otherwise, 
allowed for R and P branches; forbidden for Q branch, 
allowed for Q branch; forbidden for R and P branches, 



0, and 



are named parallel transitions (with the change in the dipole moment parallel to 
the molecular axis) and those with A£ = ±1 are named perpendicular transitions (the dipole moment changing in a 
perpendicular direction to the molecular axis). 

The three most abundant isotopologues of acetylene are 12 C 2 H 2 (C 2 H 2 ), H 13 C 12 CH (H 13 CCH) and 13 C 2 H 2 . Acety- 
lene is a linear molecule with five fundamental vibrational modes (see Figure 10). The bending modes v± and v§ are 
doubly-degenerate to a first approximation with energies of 612.871 and 730.332 cm -1 respectively (Herman et al. 
2003). The mode is infrared inactive for C 2 H 2 and 13 C 2 H 2 , but active for H 13 CCH. The v\, z^ 2 , and vibra- 
tional modes correspond to the different stretching modes of acetylene and have energies of 3372.849, 1974.316, and 
3294.839 cm -1 respectively for C 2 H 2 , and differ slightly for the isotopologues (Herman et al. 2003). Both C 2 H 2 and 
13 C 2 H 2 are symmetric and show a degeneracy, g s , in the rotational levels due to nuclear spin statistics. The para- 
levels (J = even) have g s (C 2 H 2 ) = 1 and .g s ( 13 C 2 H 2 ) = 6, while the ortho- levels (J = odd) have g s (C 2 H 2 ) = 3 and 
5«( 13 C 2 H 2 ) = 10. These two isotopologues have a null permanent dipole moment due to symmetry with respect to the 
molecular midplane. Therefore, they do not have purely rotational transitions. However, H 13 CCH is not symmetric 
with respect the molecular midplane and has a very small dipole moment along the molecular axis. The bands of 
C 2 H 2 and H 13 CCH identified are shown in the Table 2. 

HCN (and H 13 CN) has a permanent dipole moment along the molecular axis allowing purely rotational radiative tran- 
sitions. HCN has three vibrational modes: v 2l is a doubly-degenerated bending mode with an energy of 713.461 cm -1 , 
and the other two, v\ and z^ 3 , are stretching modes with energies of 3311.480 and 2096.846 cm" 1 , respectively (Maki 
et al. 2000). The vibrational energy pattern for HCN and H 13 CN is shown in Figure 15. The detected bands are given 
in Table 2. 

LINE INTENSITIES 
The line intensity as a function of temperature is given by (Jacquemart et al. 2001): 

1 87T 3 . f ow ^ 1 



S(T) 



4ire 3/ic Z(T) g, 



\R\ 2 L{J,l)e- hcE ^ /kBT 



1 



hcvo/ksT 



(Bl) 



where the dipole moment squared, \R\ is expressed in Debye 2 , v$ is the centre of the line in cm 1 , E\ ow is the energy of 
the lower level in cm -1 , <jf ow is the spin degeneracy of the lower level, L(J, £) is the rotational line strength (see below) 
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TABLE Bl 

DlPOLE MOMENTS FOR C2H2 AND HON VIBRATIONAL TRANSITIONS 





C2H2 




HCN 


Transition 


9v |R | 2 (xlO- 2 D 2 ) 


Transition 


g v |i? | 2 (xlO- 2 D 2 ) 



M*u)~ G.S.(a+) 


1 


9.81±0.11 a 


u 2 (tt) « G.S.(<r+) 


1 


3.8925 ± 0.0016 b 


U 4 + 1/5 (<T+) «-» 1/4 (7Tg) 


1 


4.904 ± 0.063 a 


2^2(0"+) *-* V2(n) 


1 


3.39 c 


^4 + ^s(o-u ) <-» ^4(7r 9 ) 

1/4 + u 5 (S) <-> ^ 4 (7r 9 ) 


1 


5.006 ± 0.062 a 


2vi(&) <-» 1^2 (tt) 


2 


7.14 c 


2 


9.599 ± 0.084 a 








2^5 (ft) <-> ^5(^11) 


1 


9.79±0.19 a 








2^(5) <-> I/ 5 (7r„) 


2 


18.86 ± 0.14 a 









Note. — See the text for a definition of g v and the ro-vibrational selection rules. 

a Jacquemart et al. (2001) 
b Malathy Devi ct al. (2005) 

c The vibrational dipolc moment has been obtained from fits to the available data in the 
HITRAN Database, 2004. Malathy Devi et al. (2005) established that the values of this database 
for the U2 <-> G.S. transition are lower than they should be. Hence we could expect the same 
for the other transitions. Therefore, we will use the same values for 2^2 (o"+) <-> 1*2 ("0 and 
2^(5) *-* V2(tt) rather than V2(tt) <-> G.S.(a- + ) (Malathy Devi et al. 2005) as an approximation. 

and Z(T) is the partition function. g v = 2 when the upper and lower vibrational levels involved in a ro-vibrational 
transition present £-type doubling, i.e., £ U p,£\ow 7^ 0; otherwise, g v = 1. In the case of g v = 2, the given value of |i?| 2 
considers the transitions having parities e and / at the same time. 

The line strength of a ro-vibrational transition, L(J,£) for a linear molecule is given by (Rothman et al. 1992; 
Hcrzberg I): 

(J+l + £)(J+l-£)/(J+l) R branch, A£ = 0, 

(2J+ l)£ 2 /J(J + 1) Q branch, A£ = 0, 

(J + £)(J — £)/J P branch, A£ = 0, 

L(J,£)={ (B2) 
{J + 2 + £A£)(J+ 1 + £A£)/2{J + 1) R branch, A£ = ±1, 

(J+ 1 + £A£)(J - £A£)(2J + 1)/2J(J+ 1) Q branch, A£ = ±1, 
(J — 1 — £A£)(J - £A£)/2J P branch, = ±1. 

The dipole moments used to calculate the spectrum of C 2 H 2 and the ro-vibrational transitions of H 13 CCH, have 
been obtained from the HITRAN Database, while the hot bands of H 13 CCH have been assumed to be identical to 
those of C2H2. The corresponding values for the v-i transition of HCN have been taken from Malathy Devi et al. 
(2005). The adopted dipole moment for HCN A^2 = l hot bands and for H 13 CN V2 transitions are similar to that of 
the vi transition of HCN. All dipole moments are shown in Table Bl. 

The partition function has been calculated by directly summing ro-vibrational levels under these conditions: 

1. For each vibrational state, at least all the rotational levels with J < J max = 65 are summed, 

2. the relative contribution to the rotational partition function of the last rotational level of each vibrational state 
must be less than 10~ 3 , which may require a sum over rotational levels having J > J max = 65, 



3. the highest considered vibrational state, v u 
in the considered transition, 



must be higher than that of the upper ro-vibrational level involved 



4. the relative contribution to the molecular partition function of u max (including the rotational partition function) 
must be less than 10~ 3 implying the possibility of summing vibrational levels with v > v max . 

For the high rotational temperatures prevailing near the stellar photosphere, condition 2 implies that the highest 
rotational level needed to be included for C2H2 in the models is between J = 116 — 130. Conditions 3 and 4 lead to 
ii m „ = 110. On the other hand, for low and intermediate Tk, the highest rotational level included to fulfill condition 
2 is well below J max = 65, e.g., for T rot ~ 400 K the error on the partition function will be 10 -3 for J ~ 32, although 



we have to sum up to J n 



65 to fulfill condition 1, and v„ 



9 (conditions 3 and 4). 
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